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ABSTRACT 

We present results from the Keck Baryonic Structure Survey (KBSS), a unique spectroscopic survey of the 
distant universe designed to explore the details of the connection between galaxies and intergalactic baryons 
within the same survey volumes, focusing particularly on scales from ^ 50 kpc to a few Mpc. The KBSS is 
optimized for the redshift range z ^ 2-3, combining S/N^ 100 Keck/HIRES spectra of 15 of the brightest 
QSOs in the sky at z ~ 2.5-2.9 with very densely sampled galaxy redshift surveys within a few arcmin of each 
QSO sightline. In this paper, we present quantitative results on the distribution, column density, kinematics, 
and absorber line widths of neutral hydrogen (H I) surrounding a subset of 886 KBSS star-forming galaxies 
with 2.0 ^ z ^ 2.8 and with projected distances < 3 physical Mpc from a QSO sightline. Using Voigt profile 
decompositions of the full Lya forest region of all 15 QSO spectra, we compiled a catalog of ^ 6000 individual 
absorbers in the redshift range of interest, with 12 <log(A^Hi)< 21. These are used to measure H I absorption 
statistics near the redshifts of foreground galaxies as a function of projected galactocentric distance from the 
QSO sightline and for randomly chosen locations in the intergalactic medium (IGM) within the survey volume. 
We find that A^hi and the multiplicity of velocity-associated H I components increase rapidly with decreasing 
galactocentric impact parameter and as the systemic redshift of the galaxy is approached. The strongest H I 
absorbers within ~ 100 physical kpc of galaxies have A^hi ^ 3 orders of magnitude higher than those near 
random locations in the IGM. The circumgalactic zone of most significantly enhanced H I absorption is found 
within transverse distances of ^ 300 kpc and within ±300 km s"' of galaxy systemic redshifts. Taking this 
region as the defining bounds of the circumgalactic medium (CGM), nearly half of absorbers with log(A^Hi) 
> 15.5 are found within the CGM of galaxies meeting our photometric selection criteria, while their CGM 
occupy only 1.5% of the cosmic volume. The spatial covering fraction, multiplicity of absorption components, 
and characteristic A^hi remain significantly elevated to transverse distances of ^2 physical Mpc from galaxies 
in our sample. Absorbers with A^hi > lO'"*'^ cm~^ are tightly correlated with the positions of galaxies, while 
absorbers with lower A^hi are correlated with galaxy positions only on ^ Mpc scales. Redshift anisotropics 
on these larger scales indicate coherent infall toward galaxy locations, while on scales of ^ 100 physical kpc 
peculiar velocities of Av ~ ±260 km s"' with respect to the galaxies are indicated. The median Doppler widths 
of individual absorbers within 1-3 of galaxies are larger by ~ 50% than randomly chosen absorbers of the 
same A^hi, suggesting higher gas temperatures and/or increased turbulence likely caused by some combination 
of accretion shocks and galactic winds around galaxies with Mhaio — 10'^ M0 at z ^ 2-3. 
Subject headings: cosmology: observations — galaxies: high-redshift — galaxies: evolution — galaxies: 
formation — intergalactic medium — quasars: absorption lines 
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1. INTRODUCTION 

Hydrogen, comprising three quarters of the baryonic mass 
of the universe, is the principal component of all luminous ob- 
jects in the universe. It is the fuel source for stars and there- 
fore for star formation. Thus, in order to understand the for- 
mation and evolution of galaxies, one must understand and be 
able to trace the inflow and outflow of this fuel. 

There exist very poor observational constraints on the 
movement of baryons in and out of galaxies. At high 
redshift in star-forming systems, it has been argued that 
the outflow rate must be similar to the star-formation rate 
dPettini et al.ll2000.) and that the inflow rate must be similar to 
the combined star-form ation rate and outflow rate ( Erb 200^ 
'Finlat or & Davll2008h . 

Recently there has been a flurry in the theoretical lit- 
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erature predicting the prevalence of accretion of cold gas 
(log(r) ^ 4.5 - 5.5 K) onto high-z galaxies via filamen- 
tary "cold flows" ("Birnboim & Dekel"2003';'Keres et al."2005'; 
Ocvirk et al. 2008; Brooks et al. 2009; Faucher-Giguere et al. 
201 It Ivan de Voort et al.ll201 laiibl) . In this model, the baryons 
stream into galaxies along the filamentary structure of the 
cosmic web, accreting onto galaxies without experiencing 
virial shocks. A wide range of predictions has been made 
concerning the efficiency of the transport of this material 
into galaxy halos, as well as its role in fueling ongoing star 
formation (van de Voort et al. 2011b; Faucher-Giguere et al. 
Ilonh . Further, there may be substantial suppression of the 
cold accretion rate caused by galaxy-scale mass outflows, ev- 
idence for which is commo nly observed in the spectra of 
high- z star-forming systems (Pettini et al.'200 lJ;IShaplev et al.l 
I2003t Adelberger et al. 2005a; Steidel et all |2010t- Map^ 
ping the gas distribution surrounding galaxies is crit ical to 
constraining these rn odels (Faucher-Giguere & Keres' '201 Ji 
iFumagalli et al.ll2"01 ll) . and would be a significant step toward 
understanding and quantifying the exchange of baryons be- 
tween the sites of galaxy formation and the nearby intergalac- 
tic medium (IGM). 

There has been a large amount of recent theoreti- 
cal examination of the nature of IGM absorbers and 
their relation to galaxies using simulations. Lya is be- 
lieved to broadly trace the filamentary large-scale structure 
jCen et al.' '1994; 'Zhang et al .1 119951: iM iralda-Escu de et"aL 
1996; Hernauist et al. 1996; R auch et a l. 1997; Theun s et al. 
1998; Dave et al. 1999; Schaye 2001) although there are in- 
dications (B arnes et al. 2011) that winds could blow spatially 
extended halos of g as which may have r ecently been observed 
both in absorption (ISteidel et al.ll2Ol"0l) and in Lya emission 
dSteidel et al.ll201 lb . There seems to be general agreement 
that galactic winds are responsible for metal absorbers in 
the IGM. Booth et al. (2010) suggest that mostly low mass 
(Mdm ^ 10'° Mf7^) galaxies must be responsible for the pol- 
lution, while IWiersma et al.l (I201Q) suggest only half of the 
metals would origina te from galaxies with Mdm ^ lO" Mq. 
IWiersma et al.l (120101) also studied the history of the ejection 
of these metals and found that half of the metals observed at 
redshift 2 were ejected during the time b etween 2 < z < 3 - 
Using cosmological "zoom-in" simulations lShen et al.l (1201 ll) 
found a "Lyman Break"-ty pe galaxy could distribute metals 
to 3 virial radii by z = 3. 'S imcod (1201 ll) recently consid- 
ered this problem observationally, finding that indeed 50% of 
the metals observed in the IG M at z ^ 2.4 were placed there 
sinc e z ~ 4.3, i.e. in 1.3 Gyr. lOppenheimer & Davd (l2008l) 
and lOppenheimer et alj (120101) studied the fate of winds us- 
ing cosmological simulations and found that while galactic 
winds are likely responsible for the metallic species seen in 
the IGM, much of the outflowing gas may be bound to galax- 
ies and may fall back in. In their simulations the recycling 
timescale scaled inversely with mass because winds emanat- 
ing from more massive galaxies experienced greater hydro- 
dynamic drag due to the increased abundance of dense IGM 
surrounding them. Further, in these simulations the largest 
source of gaseous fuel for star formation after z ^ 1 was recy- 
cled wind material. 

To date, systematic attempts to jointly study high-z 
gal axies and their i nterga l actic e nvirons ha ve been made 
bv lAdelbergeretalJ (120031 |2005a) (see also ICrighton et aP 
1201 ll) . These studies focused primarily on H I and C IV ab- 
sorption in the spectra of background QSOs whose sightlines 



probed regions covered by 2 < z < 4 Lyman break galaxy 
(LBG) surveys. This work allowed for a first glimpse of 
the distribution of diffuse gas surrounding forming galaxies 
at high redshift, and, perhaps more tantalizing, evidence for 
interactions between the IGM and galaxies during the epoch 
when galaxies a r e exp ected to be most active. Gen erally, 
lAdelbergeretalJ (120031) and lAdelberger etalj (l2005al) found 
excess H I absorption out to « 5 - 6h~'comoving Mpc (cMpc) 
of galaxies [^2 physical Mpc (pMpc) at (z) = 3.3 using the 
cosmology adopted in this paper]. C IV systems were found to 
correlate strongly with the positions of galaxies suggesting a 
causal connection. Unfortunately, these papers could not con- 
sider physical properties of the gas such as its column density 
or temperature because the data were not of sufficiently high 
quality to perform Voigt profile analysis. As such these pa- 
pers focused on the transmitted flux which could be applied 
to a wider range of data qualities. 

In this work, we provide high-accuracy analysis of the spec- 
tral regions surrounding 886 high-z star-forming galaxies as 
seen in absorption against the spectra of background hyper- 
luminous QSOs using data drawn from the Keck Baryonic 
Structure Survey (KBSS; Steidel et al, in prep). The KBSS 
was specifically designed to allow for the observation of gas 
absorption features surrounding high-redshift star-forming 
galaxies, providing unique insight into the IGM/galaxy in- 
terface at high redshift. The size and quality of the KBSS 
sample allow us to map the distribution and properties of gas 
near to individual star-forming galaxies with direct physical 
parameters such as the column density as opposed to prox- 
ies such as the equivalent width. This paper is the first in 
a series designed to study the physical properties of star- 
forming galaxies at high redshift using Voigt profile analysis 
of this data sample. The complementary analysis presented by 
iRakic et al] (1201 lb) describes the pixel statistics of the QSO 
spectra from the KBSS and the correlation of H I optical depth 
with the positions and redshifts of galaxies. 

In ^we discuss the galaxy and QSO data and present the 
Voigt profile analysis of the QSO spectra in ^ The distribu- 
tions of H I absorbers as a function of velocity, impact param- 
eter, and 3D distance are presented in ^ In ^ we consider 
the geometric distribution of the gas using the covering frac- 
tion and incidence of absorbers, ^focuses on 2D "maps" of 
the absorber distribution around galaxies. In ^we analyze 
the velocity widths of individual absorbers and their correla- 
tion with the proximity of galaxies. We discuss the results and 
their possible interpretation in ^with a brief summary of the 
paper and our conclusions in ^ 

Throughout this paper we assume a A-CDM cosmology 
with Ho = 70 km s"' Mpc"', fl^ = 0.3, and = 0.7. All dis- 
tances are expressed in physical (proper) units unless stated 
otherwise. We use the abbreviation pkpc and pMpc to indi- 
cate physical units, and ckpc or cMpc for co-moving units. 
At the mean redshift of the galaxy sample ((z) = 2.3), 300 
pkpc is 210h"'pkpc (physical) or ~ 700h"'ckpc; the age of 
the universe is 2.9 Gyr, the look-back time is 10.9 Gyr, and 
8.2 pkpc subtends one arcsecond on the sky. 

2. OBSERVATIONS 

The data presented in this paper are drawn from the Keck 
Baryonic Structure Survey (KBSS) and include a large sample 
of rest-UV (2188) and rest-optical (112) spectra of UV-color 
selected star-forming galaxies at (z) ^ 2.3. These galaxies 
were photometrically selected to lie in the foreground of one 
of 15 hyper-luminous QSOs in the redshift range 2.5 < zqso < 
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Figure 1. The redshift distribution of the 2188 galaxies in the full KBSS 
sample (dark shaded histogram) and of the subset of 886 galaxies used for 
the analysis in this paper (light shaded histogram), which are those with red- 
shifts high enough that the con'esponding wavelength of the Ly/3 transition 
is observed in the relevant QSO spectrum, and with z < zqso ~ 3000 km s"' . 
The dark vertical line segments mark the redshifts of the 15 KBSS QSOs. 

2.85 for which we have obtained high-resolution, high signal- 
to-noise ratio (S/N) echelle spectra. 

The redshift range of this survey has important sig- 
nificance in the history of the universe - it coincides 
with the p eak of both universal star formation (see 
iReddv et al.l l2008h and super-massive black hole growth 
(see iRichards et al] I2006h . Spectroscopic observations of 
star-forming galaxies during this epoch commonly exhibit 
signatures of strong outflowing wind s (Petti ni et al . 2001; 



Shaplev et al] l2003t lAdelberger et al.1 l2005al: ISteidel et d] 
2010) . At the same time, the baryonic accretion rate 



onto galaxies is predic ted to be near its peak at z ^ 2.5 
(e.g.. Van de Voort et alj|201 Ibt iFaucher-Giguere et alj|20Tlt 
lyan de Voort et al. 20 1 lai ) TTius. the signatures of galaxy for- 
mation within the IGM should be at their peak as well. 

The redshift range 2.0 ^ z ^ 2.8 offers a number of prac- 
tical advantages as well: first, the rapidly-evolving Lya for- 
est has thinned enough to allow measurements of individual 
H I systems and to enable the detection of important metallic 
transitions falling in the same range (notably, O VI AA1031, 
1036); second, the astrophysics-rich rest-frame far-UV be- 
comes accessible to large ground-based telescopes equipped 
with state-of-the-art optical spectrographs; third, the rest- 
frame optical spectrum is shifted into the atmospheric trans- 
mission windows in the near infrared, allowing for obser- 
vations of a suite of diagnostic emission lines (Ha, H/3, 
[O II] AA3726,3729, [O III]AAA4363,4959,5007, and [S II] 
AA67 17,6732) neatly packed into the J, H, and K bands. 

2.1. The Galaxy Sample 

The KBSS galaxy sample, a subset of which is used in this 
paper (Figure[T]i, is described in detail by Steidel et al (2012); 
here we present a brief summary. Galaxies were selected for 
spectroscopy using their rest-frame UV colors ( i.e., LBGs) 
according to the criteria outlined by 'Steidel et al.' (2004) and 
lAdelberger et al.. (.2004.) for z ~ 2 and bv . Steidel et al..a003l) 
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Figure 2. The number of galaxies as a function of physical impact parameter 
Aran for the sample which has appropriate redshifts to be used in this work 
(Ught shaded histogram in Figure[T] see Table[T). The decline in the number 
of galaxies at Dtian ~ 2 pM pc is d ue to the typical survey geometry of the 
KBSS fields as described in PTT] 

for z ^ 3. In combination, these criteria have been devised to 
efficiently select star-forming galaxies over the redshift range 



1.5, 



, 3.5. In total, the KBSS galaxy sample includes R : 



1000 (^5-lOA) rest-UV spectra of 2188 star-forming galaxies 
obtained at the W.M. Keck Observatory using the Keck 1 10m 
telescope and the blue arm of the Low Resolution Imagi ng 
Spectrometer (LRIS-B: IOkeetan[T99llSteidel et al.ll2004l ) 

Most of the observations were conducted using a 
400 lines mm"' grism (blazed at 3400 A in first order) in com- 
bination with a dichroic beamsplitter sending all light short- 
ward of ^ 6800 A into the blue channel, where the wavelength 
coverage for a typical slit location was ~ 3 1 00 - 6000 A with 
a resolving power R ~ 800 using 1'.'2 wide slits ( ISteidel et alj 
12010). Some of the spectra were obtained using the d560 
beamsplitter (beam divided near 5600 A) together with a 
600 lines mm"' grism (4000 A blaze), typically covering 
~ 3400 - 5600 A with a resolving power 7? ~ 1 300. Observa- 
tions obtained after July 2007 made use of the Cassegrain At- 
mospheric Dispersion Corrector, thus minimizing the effects 
of differential atmospheric refraction over the spectral range 
of interest. Wavelength calibration was accomplished using 
Hg, Cd, Zn, and Ne lamps, with zero point corrections based 
on night sky emission lines on each individual exposure. '"* 

The exposure times allocated to individual galaxies ranged 
from 1.5-7.5 hours depending on the number of separate 
masks containing the same target; typically galaxies were ob- 
served on either 1 or 2 masks, each mask receiving a total 
integration of 1.5 hours. Further details on the selection, ob- 
serving strategy, an d data reduction fo r KBSS galaxies are 
presented elsewhere (ISteidel et alj|2010l Steidel et al. 2012). 

The rest-UV spectra of LBGs are dominated by the con- 
tinuum emission of O and B stars, over which are super- 

'** An extensive discussion of the LRIS/KBSS wavelength cahbration 
errors (including slit illumination and atmospheric dispersion) is given in 
[Steidel et al.. aOlQ) . 
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Table 1 

KBSS Central QSOs and Foreground Galaxy Samples 



Name 


RA 


Dec 


ZQSO'' 


\ b 
'»min 


Zgai range 




S/N Lyo'' 


S/NLy/3 


QOlOO+130 (PHL957) 


01:03:11.27 


+13:16:18.2 


2.721 


3133 


2.0617- 2.6838 


47 


77 


50 


HS0105+1619 


01:08:06.4 


+16:35:50.0 


2.652 


3230 


2.1561- 2.6153 


53 


127 


89 


Q0142-09 (UM673a) 


01:45:16.6 


-09:45:17.0 


2.743 


3097 


2.0260- 2.7060 


65 


71 


45 


Q0207-003 (UM402) 


02:09:50.71 


-00:05:06.5 


2.872 


3227 


2.1532-2.8339 


46 


82 


55 


Q0449-1645 


04:52:14.3 


-16:40:16.2 


2.684 


3151 


2.0792- 2.6470 


50 


73 


41 


Q0821+3107 


08:21:07.62 


+31:07:51.17 


2.616 


3239 


2.1650-2.5794 


37 


50 


33 


Q1009+29 (CSO 38) 


10:11:55.60 


+29:41:41.7 


2.652= 


3186 


2.1132-2.6031 


36 


99 


58 


SBS1217+499 


12:19:30.85 


+49:40:51.2 


2.704 


3098 


2.0273- 2.6669 


43 


68 


38 


HS 1442+2931 


14:44:53.67 


+29:19:05.6 


2.660 


3152 


2.0798- 2.6237 


46 


99 


47 


HS1549+1919 


15:51:52.5 


+19:11:04.3 


2.843 


3165 


2.0926- 2.8048 


54 


173 


74 


HS 1603+3820 


16:04:55.38 


+38:12:01.8 


2.551' 


3181 


2.1087-2.5066 


37 


108 


58 


Q 1623+268 (KP77) 


16:25:48.83 


+26:46:58.8 


2.5353 


3126 


2.0544- 2.4999 


133S= 


48 


28 


HS 1700+64 


17:01:00.6 


+64:12:09.4 


2.751 


3138 


2.0668-2.7138 


im 


98 


42 


Q2206-199 


22:08:52.1 


-19:43:59.7 


2.573 


3084 


2.0133- 2.5373 


45 


88 


46 


Q2343+125 


23:46:28.30 


+12:48:57.8 


2.5730 


3160 


2.0884- 2.5373 


84*5 


71 


45 



" The redshift of the QSO 
The minimum wavelength covered in the HIRES QSO spectrum 

The number of galaxies in our LRIS survey with Aran < 3 pMpc and spectroscopic redshifts in the correct range given in the 
previous column. 

The average signal to noise ratio per pixel of the QSO spectrum in the wavelength range pertaining to COM Lya and Ly/3 
absorption. 

The redshift of this QSO was revised after the fitting of the HIRES spectrum was completed. The redshift assumed for the 
Av = -3000 km s"' is 2.6395 

' The redshift of this QSO was revised after the fitting of the HIRES spectrum was completed. The redshift assumed for the 
Av = -3000 km s"' is 2.5420 

8 The photometry and spectroscopy in fields Q1623, HS1700, and Q2343 cover larger areas than a single LRIS footprint. As 
such, they represent the fields in which we sample the galaxy-IGM connection at Aran > 2 pMpc. 



posed numerous resonance absorption lines of metallic ions 
and H I. The H I Lya line at 1215.67 A may be seen in emis- 
sion or absorption (and often in both). The absorption fea- 
tures arise in cool interstellar gas in the foreground of the OB 
stars; they are most commonly observed to be blue-shifted 
by 100 - 800 km s"' with respect to the systemic velocity of 
the stars, as measured from either rest-frame optical nebular 
emiss ion lines or stellar photospheric lines i n stacked spectra 
(see Pettini et al."2001'; 'S haplev et aLl l2003^, 'Ade lberger et alj 
|2005a; Steidel et al. 2010) or from the redshift-s pace symme - 
try of Lya absorption in the nearby IGM (iRakic et al.ll201 lah . 
Common lines observed include: O VI AA1031,1036, Si II 
A 1260, A 1526, Si IH-O I A1303 (blend),Si III A 1206, Si IV 
AA1393,1402, N V AA1238,1242, C II A1334, C III A977, and 
C IV AA1548,1550. The profile of the Lya emission or ab- 
sorption line is modulated by the optical depth of the material 
closest to the systemic velocity of the stars, which has been 
shown t o correlate most significantly with the baryonic mass 
(ISteidel et al. 2010) and the physical size (Law et al., in prep) 
of the galaxy. 

The sample used in this study includes 886 galaxies within 
the redshift range where at minimum Lya and Ly/3 are ob- 
served in the HIRES spectrum and with redshifts placing them 
at least 3000 km s"' blue-ward of the redshift of the QSO. The 
latter criterion was selected to avoid proximate systems that 
originate within material ejected from the QSO itself and/or 
the region affected by its ionizing radiation field. 

A typical galaxy in the spectroscopic survey has a bolo- 
metric luminosity of 2.5 x 10" Lq jReddv et all I2008L 
201 1), a stai--formationrate (SFR) of 30 Mf7> yr"' (lErb et alJ 
2006bl a stellar age of - 0.7 Gyr (Erb et al."2006c^, and a 
gas-phase metallicity of ~ 0.5 Zq (Erb et al. 2006a). The 
gala xies inhabit dark matter halos of avera ge mass ^ 10'^ 
M0 (lAdelbergeretal.ll2005bt IConrov et al.ll2008> Trainor & 
Steidel, in prep; Rakic et al, in prep) and average dynamical 



masses of ^ 7 x 10'" M0 (Erb et al."2006c) and generally ex- 
hibit dispersion-dominated kinematics (Law et al. 2009). The 
luminous parts of the galaxies are do minated by baryons , 
typical ly half st ars and half cold gas dShaplev et alJ 120051; 
Erb et al]l2006cl) . with half-light radii of ~ 2 pkpc jLaw et all 
201 1). The spectroscopic sample includes objects with appar- 
ent magnitudes TZ < 25.5, where TZ is equivalent to mAB(6830 
A). At the mean redshift of the sample ((z) = 2.30) the faint 
limit corresponds to a galaxy of 0.25Luy (iReddv & Steidell 
I2OO9I) . 

The redshift distribution of the galaxy and QSO sample is 
presented in Figure[Tl and the distribution of physical impact 
parameters, Dtran, between the galaxies and the QSO lines of 
sight is shown in Figure |2] 

At present, 1 12 galaxies in the full KBSS sample have been 
observed spectroscopically in the near-IR using using NIR- 
SPEC (McLean et al. 1998) on the Keck II telescope. 87 of 
these galaxies lie in our chosen redshift interval. The NIR- 
SPEC target selec tion, data, and reductions are discussed in 
lErb et all (l2006blB . The NIR redshifts, generally based on 
the Ha emission line, are estimated to be accurate to ^ 60 
km s-i or cr, ~ 0.0007 at z - 2.3. 



2.2. Measured and Calibrated Redshifts 

Because the most prominent features in the UV spectra of 
star forming galaxies are not at rest with respect to a galaxy's 
systemic redshift, Zgai, corrections must be applied to avoid 
substantial systematic redshift errors. The velocity peak and 
centroid of the Lya emission line, when present, tend to be 
redshifted with respect to Zgai by several hundred km s~', 
while the strong UV absorption features (zjs) tend to be sim- 
ilarly blue-shifte d with respect to z„at dShaplev et alj 120031 ; 
lAdelbergereTal] 120031; ISteidel et al.ll2010l) . These observa- 
tions are generally interpreted as strong evidence for the pres- 
ence of galaxy-scale outflows. 
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Here we adopt estimates of galaxy systemic r edshifts 

S Zgai) co mputed in the manner pro posed by [A delberg er et aTl 
2005ab an d later updated by Steideletal. (2010^ and 
Rakic et alJ ([2011 ah . [Adelberger et al. (2005a) and 
Steidel et al.l (|2010^ analyzed the subset of the UV sam- 



ple for which both rest-UV and rest-optical spectra had been 
obtained. They measured the average offset between redshifts 
defined by Ha emission versus Z]_ya and z\s to estimate 
average corrections. The Ha line traces the ionized gas in 
star-forming regions and is therefore a reasonable proxy for 
the systemic velocity of the stars, which are more difficult 
to measure due t o the weakn e ss of t he UV photospheric 
absorption lines. iRakic et aTl (1201 lal) used the QSO and 
galaxy data set presented here and calibrated velocity offsets 
appropriate for various classes of LBGs by insisting that the 
average IGM Lya absorption profiles should be symmetric 
with respect to galaxy redshifts. In both cases, the offsets 
represent those of the ensemble while in reality there is 
some scatter between individual objects even if their spectral 
morphology is similar However, as we will demonstrate, 
the adopted Zgai must be generally quite accurate in order to 
produce the trends described below. 

The formulae used for estimating Zgai from ZLya and zis 
measurements are reproduced below. For galaxies with Ha- 
based redshifts (87/886), we set Zgai = zhq- For galaxies which 
have measured zis with or without the presence of Lya emis- 
sion (691/886), 



where 



_ Avis 

Zgal.IS — ZlS+ ^ (1 +Zis) , 



Avis = 160 km s ' 



(1) 
(2) 



is the velocity shift needed to transform the observed red- 
shift into its systemic value, zis is the measured redshift from 
the centroids of interstellar absorption lines, and Zgai.is corre- 
sponds to the estimated systemic redshift of the galaxy. 

For galaxies which have redshifts measured only from Lya 
in emission (90/886), we compute the redshift as 



^gal,Lya — ^Lya 



■(l+ZLyo), 



where 



A 



Vhya '■ 



-300 km s" 



(3) 



(4) 



is the velocity shift needed to transform the observed redshift 
into its systemic value, Zhya is the measured redshift from 
Lya, and Zgai.Lya is adopted systemic redshift. 

For galaxies with measurements of both zis and ZLya, we 
verify that zis < Zgai.is < Zhya- If the corrected absorption red- 
shift is not bracketed by the two measured redshifts (18/886 
galaxies)'^, then we use the average of zis and ZLya^ 



Zgai : 



ZLya+ZlS 



(5) 



The residual redshift errors have a significant impact on our 
ability to interpret the kinematic information in the data; thus, 
their amplitude will be important to consider in the exam- 
ination of the line-of-sight distribution of H I. ISteidel et al.l 

In these cases, zis and Z]_ya have very similai' values and therefore taking 
their average generally provides a systemic redshift with less potential error 
than in those galaxies for which an average shift (as in equations [T]and[3) is 
required. 



( 120101) found this method generally corrects the redshifts to 
within ~ 125 km s"' of the systemic velocity. 

2.3. QSO Observations 

The 15 hyper-luminous {my ~ 15.5- 17) QSOs in the cen- 
ter of the KBSS fields (Table[Tll were obser ved with the High 
Resolution Echelle Spectrometer (HIRES: IVogt et al.|[T994h 
on the Keck I telescope. All available archival data for these 
15 QSOs have been incorporated, including data taken with 
UVES (Dekkeret al. 2000) on the VLT for Q2206-199 and 
Q2343-I-125. We obtained additional HIRES observations in 
order to reach a uniformly high S/N ratio over the spectral 
range of primary interest, 3100-6000 A. The final HIRES 
spectra have R ~ 45,000 (FWHM- 7 km s"'), S/N - 50-200 
per pixel, covering at least the wavelength range 3100 - 
6000A with no spectral gaps. The significant improvement 
in the UV/blue sensitivity of HIRES resulting from a detec- 
tor upgrade in 2004 enabled us to observe Ly/3 A 1025. 7 down 
to at least z = 2.2 in all 15 KBSS sightlines, and to signifi- 
cantly lower redshifts in many (see Table 1.) The additional 
constraints provided by Ly/3 (and in many cases, additional 
Lyman series transitions) allow for much more accurate mea- 
surements of H I for Nm= 10'"*- 10'^ cm"^ (see ®; this is 
particularly important since these column densities are typi- 
cal of H I gas in the CGM at these redshifts (see 

The QSO spectra were reduced using T. Barlow's MA- 
KEE package which is specifically tailored to the reduction 
of HIRES data. The output from MAKEE is a wavelength- 
calibrated'^ extracted spectrum of each echelle order, cor- 
rected for the echelle blaze function and transformed to vac- 
uum, heliocentric wavelengths. The spectra were continuum- 
normalized in each spectral order using low-order spline in- 
terpolation, after which the normalized 2-D spectrograms 
were optimally combined into a single one-dimensional, 
continuum-normalized spectrum resampled at 2.8 km s"' per 
pixel. 

Three of the KBSS QSO spectra contain a damped Lyman 
a (DLA) system (A^hi> lO^" -' cm"^), and three contain a sub- 
DLA (10''' ^ A^Hi< 10^"^ cm"^). Special care must be taken 
in the continuum fit to the regions surrounding these systems, 
as the damping wings of the absorption lines extend for thou- 
sands of km s~' from line center In these regions, we care- 
fully fit a Voigt profile to the core of the absorber and adjusted 
the original continuum fit so that the Voigt profile produced a 
good fit to the damping wings, as shown in the top panel of 
Figure [3] The final Voigt profile fit is divided into the true 
spectrum, resulting in a new spectrum where the wings of the 
DLA have been removed, as illustrated in the bottom panel 
of Figure [3] The re-normalized spectrum can then be used to 
fit additional absorption systems superposed on the damping 
wings. 

The redshifts of the QSOs are measured from rest-frame op- 
tical emission lines using lower-dispersion NIR spectra. The 
details of this procedure and the expected errors in the QSO 
redshifts are reported upon in Trainor & Steidel (2011, sub- 
mitted). The precise QSO redshifts do not affect our analysis. 

3. ANALYSIS OF QSO ABSORPTION SPECTRA 

The wavelength calibration of the HIRES spectra introduce negligible 
eiTor into our analysis. The HIRES spectra have calibration errors less that 
0.5 km s"'. 
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Figure 3. A demonstration of our treatment of the continuum surround- 
ing damped-LycK systems. Top: In black, the continuum-normaHzed HIRES 
spectrum of Q2343+125 showing ±7000 km s"' surrounding the DLA. The 
(red) dashed line corresponds to the Voigt profile of the DLA centered at 
z = 2.4312 with log(A'Hi) = 20.4. Shown in the light (blue) curve is the error 
spectrum. Bottom: The HIRES spectrum of the same QSO with the DLA 
profile divided out. The new error spectrum accounting for the DLA profile 
division is shown by the light (blue) curve. The new continuum (with the 
DLA divided out) is shown in the dashed (red) curve. 

The process of accurately measuring H I in the Lya forest 
of QSO spectra is compHcated by the saturation of moderately 
strong absorbers and the blending of H I features with other 
H I or with lines of metallic species that happen to fall in the 
forest region. 

Our analysis includes a full Voigt-profile decomposition of 
the Lya forest from the lowest redshift for which Ly/3 is avail- 
able in each spectrum up to 3000 km s"' blue-ward of the 
QSO redshift; Table 1 shows the relevant redshift range used 
for each QSO in the sample. The cut-off at the high redshift 
end is to avoid H I systems which could be ejected from and/or 
ionized by the QSO itself; the low-redshift cut is necessary 
due to the high frequency of systems with A^hi ^ lO'**'^ cm"^ 
that will be saturated in Lya'^. For saturated systems, the 
fit to Lya is degenerate between an increase in column den- 

Note that the exact saturation point for the Lyo line and for any other 
line depends on the line width (i.e. the Doppler parameter) and thus the 
internal properties of the absorbing cloud. 



sity or an increase in the width of the line, b^. The best way 
to resolve this degeneracy is to measure higher-order Lyman 
lines where decreasing oscillator strengths allow accurate A^hi 
determination. 

Simultaneous fits were made to as many Lyman lines as 
were both (1) available in the observed spectral range and (2) 
needed to measure an unsaturated and uncontaminated profile 
in the highest-order line. The exact number of higher-order 
Lyman lines used therefore depended upon both the redshift 
of the absorber and the degree of contamination in the spec- 
tral region containing the higher-order Lyman series absorp- 
tion features. Higher-order lines, whether saturated or unsatu- 
rated, were used whenever doing so provided additional con- 
straints on the overall fit. 

Example Voigt profile fits to the H I absorption in regions 
surrounding the systemic redshift of five galaxies from our 
sample are shown in Figure |4] Note that for galaxies with 
redshifts significantly larger than z ^ 2.2, many higher-order 
Lyman transitions can be measured. 

To facilitate the fitting of the Lyman a forest, we devel- 
oped a semi-automatic line-fitting code. Briefly described, 
the code works with ^ 1500 km s"' sections of spectrum at 
a time'**, fitting to Lya and as many higher-order lines as are 
accessible within the HIRES spectrum at the redshift of the 
H I systems being fit. The algorithm first searches for systems 
by cross correlating a template hydrogen absorption spectrum 
(i.e., a single non-saturated H I absorption component) with 
the HIRES spectrum. Peaks in the cross correlation are taken 
as initial estimates of the centroids of absorption lines. We 
fit Gaussians to these lines to estimate column densities and 
Doppler parameters. Residual absorption features (i.e., those 
inconsistent with Lya) in the Lya portion of the spectrum are 
assumed to be metal lines. Residual absorption in the higher- 
order Lyman series sections of the spectra are assumed to be 
lower-redshift H I systems. The fits begin with Lya at the 
high-redshift end of the range so that their higher-order Ly- 
man absorption can be flagged as a known contaminant for 
fitting lower-redshift H I absorbers. 

Once estimates of the locations, column densities, and 
Doppler parameters of all the absorption lines are complete, 
they are input into the minimization code VPFIT'^ written 
by R.F Carswell and J.K. Webb. VPFIT simultaneously fits 
all transitions of H I as well as the specified contaminating 
lines. The results are checked by eye, alterations made where 
the fit is inappropriate, and the process is repeated iteratively 
until a good fit (reduced ~ 1) is achieved. At this point, the 
multiple sections of spectrum are spliced together until a full 
fit to the forest is achieved. 

It should be noted that the Voigt profile fit to a spectrum 
does not represent a unique solution. In this work, we fit each 
set of absorbers with the minimum number of components, 
adding additional components only when they significantly 
improve the x^- Median errors in A^hi reported by VPFIT are 
0.07 dex for absorbers with 13 < log(A^Hi) < 14 and 0.03 dex 
for absorbers with 14 < log(A^Hi) < 16; however in many cases 
the systematic errors will exceed these values. The largest 
source of error in our fits to low column-density systems is 

Note, the choice of 1500 km s"' was made largely out of convenience. 
This velocity window is the largest window which can be easily displayed to 
check the goodness of fit and for which the number of components allows for 
a reasonable VPFIT run time. 

" http://www.ast.cam.ac.uk/ rfc/vpfit.html; © 2007 R.F. Carswell, J.K. 
Webb 
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Figure 4. Example fits to the QSO data surrounding the redshifts of galaxies in our sample. Displayed in black are the continuum-normalized HIRES spectra 
showing the Lyman series transitions within ±700 km s"' of the systemic redshift of 5 galaxies with redshifts as indicated. Over-plotted in color are Voigt profile 
decompositions for H I absorption systems within ±700 km s"' of the galaxy redshift. Successive rows illustrate the fit to Lya, Ly/3, etc. Absorption in the 
HIRES spectra that does not appear in the colored fit corresponds to absorption from metallic species or from H I at a redshift far from that of the galaxy. Note 
that for all galaxies in our sample, the QSO spectra cover the Lyo and Ly/3 transitions near the galaxy redshifts, and for the higher-redshift galaxies, many more 
transitions in the Lyman series can be measured. 



uncertainty in the continuum level; for high column-density 
absorbers, it is the possibility of unrecognized sub-component 
structure. 

The complete decomposition of the Lya forest in these 15 
lines of sight includes Voigt profile fits to > 5900 distinct 
H I absorption systems with A^hi> 10'^ " cm"^, making it the 
largest absorber catalog ever compiled at these redshifts. It 
increases by an order of magnitude the number of intermedi- 
ate A^Hi absorbers measured with the additional constraint of 
higher-order Lyman lines. 

4. CIRCUMGALACTIC H I 

In the following sections, we discuss the statistics of indi- 
vidual H I absorption systems with respect to the redshifts and 
transverse positions of galaxies. We do not uniquely "assign" 
each absorber to a specific galaxy or vice versa. Instead we 
rely on comparisons between the absorption measured close 
to galaxies with that typical of "random" locations^" in the 
IGM. This allows us to quantify the significance of any ap- 
parent correlation with galaxies. 

In principle, it would be preferable to compare absorbers 
found close to galaxies with those found in IGM locations 

A "random" location, as described later in this section, is in effect a 
random redshift from our galaxy redshift catalog and an independent random 
position on the sky from our galaxy position catalog. 



known to be far from galaxies; however, because the galaxy 
sample is spectroscopically incomplete compared to our pho- 
tometrically selected targets,^' we do not measure the red- 
shifts of all galaxies in our survey volume. As a consequence, 
we do not know which locations in 3D space do not have a 
nearby galaxy. Thus, we can only compare locations near to 
galaxies with random locations in the IGM (irrespective of the 
positions of galaxies). 

In order to reproduce the absorber distributions for a typi- 
cal place in the IGM, we compiled a catalog of 15,000 random 
locations (in both redshift space and on the plane of the sky). 
The redshifts are drawn from the actual galaxy redshift list 
and therefore reproduce the typical IGM absorption associ- 
ated with the redshift ranges covered by our galaxy sample. 
With each redshift, we also associate a randomly drawn QSO 
and an impact parameter to the QSO sightline from the real 
galaxy impact parameter list. We can then study the distri- 
bution of absorption systems around these 15,000 random lo- 
cations and thereby understand how the presence of a galaxy 
alters the distributions. 

Below we consider the distribution of H I surrounding 
galaxies: first along the line of sight, then on the plane of the 

^' Our ph otometric sample i s also incomplete at the faintest apparent mag- 
nitudes (see lReddv et alj|2008l) . and no attempt is made to include galaxies 
with -R > 25.5. 
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Figure 5. The velocity-space distribution of H I absorption systems with respect to the systemic redshift of galaxies, normalized by the number of galaxies in the 
sample. Absorbers with logCWni) > 13 and within 1 pMpc of the sightline to a QSO are included. The solid histogram represents the distribution of H I around 
galaxies, whereas the hatched histogram represents the average absorber density near randomly-chosen redshifts drawn from our galaxy redshift distribution. 
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Figure 6. As in Figure|5] where the histogram is A'm-weighted as described in 

sky, and finally as a function of 3D distance. These measure- 
ments are used to determine the relevant velocity and trans- 
verse scales of circumgalactic H I. 

4.1. Velocity-Space Distribution ofHl Near Galaxies 

The properties of H I gas near galaxies in our sample can be 
quantified in several ways. First, we consider the line of sight 
velocity distribution of absorbers relative to the redshifts of 
galaxies. Shown in Figure|5]is the velocity distribution of all 
absorbers with A^hi> lO'-' cm"^ within 1400 km s"' in redshift 
and 1 pMpc in projected distance from a galaxy. We define the 
velocity offset, Av, of an absorber, 

^^^(Zabs-Zeal)c 
1+Zgal 

where z^bs is the absorption system redshift an d Zga i is the 
adopted systemic redshift of the galaxy from ^2.21 With 
this definition, absorbers blue-shifted with respect to galax- 
ies have negative Av. 

We define P(Av,Dtran) as the number of absorbers per 
galaxy at a given Av and within the specified range of impact 
parameters Dtran- The solid histogram in Figure |5] represents 



1.6 




|Av)| [km s '] 

text. 

the distribution of H I around galaxies, whereas the hatched 
histogram shows the average number of absorbers expected 
relative to randomly chosen redshifts, as described above. The 
number of absorbers is clearly higher near galaxies, with an 
excess peaking near Av = (the galaxy systemic redshift) and 
confined to Av ~ ±700 km s"' . 

Figure |6] shows a similar pattern; in this case, each absorp- 
tion system is weighted in proportion to log(A^Hi) such that 
high-column density systems contribute more significantly to 
the histogram. Each absorber with log(A^Hi)> 13 contributes 
[log(A^Hi) -13] to the histogram, which is normalized by the 
number of galaxies considered. Taken together, Figures|5]and 
|6] indicate that there is an increase near galaxies of both the 
number and the column density of H I absorbers. The nar- 
row peak of the Av distribution has an apparent half-width of 
-300 km s"', while the full excess extends to ~ ±700 km s 
(most clearly shown in the right-hand panel of Figure|6]l. 

As has been arg u ed by Shaplev et al. (2003); 
lAdelberger et"!!] (l2003h ; ISteidel et al.1 (120101) ; iRakic et al] 
(1201 lal) . these velocity distributions are also useful for 
checking our redshift calibration; their symmetry about 
Av= km s"' is a sensitive probe of systematic errors in our 
galaxy systemic redshift calibration, while the width of the 
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Figure 7. The median value of Nm as a function of the velocity offset of 
absorbers with respect to those galaxies with impact parameters Dam< 300 
pkpc. Asterisks represent the median value of A'hi, dark vertical bars ai'e 
the \-cr dispersion in the median determined via the bootstrap method. The 
dark horizontal line is the median value of Whi in the random distribution. 
The light shaded boxes are the bootstrapped symmetric l-cr dispersion in the 
median values of the samples drawn from the random distribution. The top x- 
axis shows the c onve rsion between velocity offset and Hubble distance. See 
Equation|2]and j|4.3l for fuilher discussion. 
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Figure 8. Nm as a function of Av for absorbers with Dtran< 100 pkpc. The 
dashed (red) line marks the position of logCWm) = 14.5. The significance of 
this column density threshold is discussed in 

distribution provides an upper limit on the random errors. 
A Gaussian fit to Figure |5] yields a mean Av = 18 ± 26 with a 
standard deviation of 308 ± 30 km s~' . 

An alternative method of quantifying the column density 
dependence in Figure|6]is to examine the distribution of A^hi 
a function of Av. Figure|7]shows the median A^hi as a function 
of Av for all absorbers within Z)tran< 300 pkpc of a galaxy-^. 

The width is an upper limit because it includes both redshift eiTors and 
any peculiar velocity of absorbers relative to galaxies. 

We choose 300 pkpc as our distance cut because the majority of the 



The asterisks indicate the median value of A^hi for each bin in 
velocity space. The dark vertical lines represent the disper- 
sion in the median computed through the bootstrap method^'*. 
The value of the median column density across all velocity 
bins for the random distribution is shown as the dark hori- 
zontal bar; the light shaded contours are the l-a bootstrapped 
dispersions in the median of the random sample, where we 
consider samples of the same size as those from the real dis- 
tribution. Note that there is an enhancement by a factor of 
~ 3 - 10 in the median A^hi out to | Av| ~ 300 km s"' rela- 
tive to galaxy redshifts. However, for | Av| > 300 km s"', the 
measurements are consistent with random places in the IGM. 

Shown in Figure |8] are the individual measurement of A^hi 
as a function of Av for all absorbers within 1000 km s"' of a 
galaxy within Z)tran< 100 pkpc of a QSO sightline. Notably, 
the higher-A^Hi absorbers cluster strongly near the galaxy sys- 
temic redshift with a full width of ±300 -400 km s"'. 

In summary, Figures|5]-[8]clearly illustrate an enhancement 
in both the column density and number of absorption systems 
near the systemic velocities of galaxies. The most significant 
enhancement of column density is seen within ± 300 km s~' 
(Figure |7]l, but there is a higher number of absorbers out to at 
least ± 700 km s"' (Figures|5]and|6]l. 

4.2. Transverse Distribution of Absorbers 

In addition to the strong velocity alignment of absorption 
systems with the systemic redshift of nearby galaxies, there 
is also a significant increase in the column densities, A^hi, of 
individual absorbers with decreasing projected (or transverse) 
distance between the galaxy and the line of sight, Dtran- Here, 
we suppose that, as an ensemble, these galaxies show simi- 
lar circumgalactic absorption signatures. Therefore, because 
galaxies fall at various discrete impact parameters from the 
QSO line of sight, we can combine the information from each 
galaxy to make a sparsely sampled map of the absorption as a 
function of Dtran relative to the ensemble galaxy. 

We introduce two related statistics designed to trace the 
change in column density as a function of Z)tran- Recalling that 
±700 km s"' encompassed the bulk of the "excess" absorp- 
tion 04.11) . for each galaxy we define max(A^Hi, 700km s"') 
to be the value of log(A^Hi) for the strongest absorber 
with |Av|<700 km s~' of the galaxy systemic redshift. 
The left panel of Figure |9] shows the median value of 
max(A^Hi,700km s"') as a function of impact parameter. A 
second statistic is the logarithm of the total Nm, sum(A^Hi), of 
all absorbers with |Av|<700km s"'. The statistics of the me- 
dian value of sum(A^Hi,700km s"' ) versus Dtran is shown in the 
right-hand panel of Figure |9l Generally, the values of these 
two statistics are quite similar because the A^hi in most veloc- 
ity windows is dominated by the single highest-A^Hi absorber. 
We consider the sum because it is most easily compared to re- 
sults of numerical simulations, as it does not require the fitting 
of Voigt profiles to simulated data and can instead be com- 
pared to a simulation "collapsed" along the line of sight. 

Figure |9] clearly demonstrates that both max(A^Hi) and 
sum(A^Hi) increase rapidly as one approaches a galaxy. In the 
bin corresponding to the smallest impact parameters, Z)tran< 

excess absorption is found within that zone, as desciibed in j|4.2l 

^'^ To bootstrap the dispersion in the median, we draw (with replacement) 
sets of data with the same size as the real distribution. For each set the me- 
dian is computed; 100 such data sets are evaluated, and then we take the l-cr 
symmetric bounds on the distribution of the medians, i.e. the 16th and 84th 
percentiles. 
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Figure 9. The log column densities of the strongest A'hi absorbers as a function of transverse distance. On the left we consider the max(A'Hi) statistic, logCWm) 
of the single strongest absorber per galaxy with | Av| < 700 km s"' . On the right is the sum(A'Hi) statistic, the log of the sum of the Whi of all the absorbers within 
I Av|< 700 km s"' . Asterisks represent the median value of the considered statistic in a given bin of Aran- Dark vertical bars are their dispersions. The horizontal 
position of the asterisks represent the median D(ran of the galaxies in that bin. The number of galaxies in each bin is indicated at the bottom of the plot. The dark 
horizontal line is the median value drawn from the random distribution. The light shaded boxes are the bootstrapped symmetric l-cr dispersion in the median 
values of the samples drawn from the random distribution. The bin size is 100 pkpc for absorbers with Z)tran< 1 pMpc and 200 pkpc for those with Z)tian> 1 
pMpc. We increase in binning to reduce the shot noise in the bins at Dtran > 2 pMpc which have fewer galaxies due to the limited field-of-view of LRIS. 
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Figure 10. Same as Figure|9]but for the maximum column density absorber 
within ±300 km s"'. Changing the velocity interval considered with the 
maxC^Hi) statistic has little effect on the observed trends. 

100 pkpc, the median value of max(A^Hi) is more than two 
orders of magnitude higher than that of a random location. 
Moving outwards, the median value decreases with increas- 
ing Aran to 300 pkpc, at which point the statistic "plateaus" 
and remains significantly higher than the random sample out 
to 2 pMpc. The plateau value in the galaxy-centric sample 
is max(A^Hi)— lO''*^ cm"^, while that of the random distribu- 
tion is max(A^Hi)— 10'^ ' cm"^. As we will discuss in ^4.2.1! 
max(A^Hi) begins to decline for Dtran> 2 pMpc. 
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Figure 11. The individual measured values of max(A'Hi, 300km s"') for 
galaxies with Dtran < 300 pkpc. 

A relevant question concerns the dependence of these statis- 
tics on the size of the velocity window considered. In Fig- 
ure |9| we considered the maximum column density absorber 
within ±700 km s"' of the systemic velocity of each galaxy. 
This corresponds to the full width of the velocity distribution 
shown in Figure |6| However, it is clear from Figures |6| and 
[Tjthat the majority of the excess strength of absorption falls 
within ±300 km s"' , especially for those systems with small 
impact parameters. FigurefTOlshows max(A^Hi,300km s"'); the 
trends are similar, though in the more restricted velocity win- 
dow the peak on small scales is higher relative to random IGM 
locations-the median value in the first bin is 3 dex higher 
than the random redshift sample, and the extended floor of 
absorption is increased to at least .5 dex above the median of 
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Figure 12. hya absorption within ±1000 km s"' of the systemic redshift of 
the 10 galaxies within Z)tran<100 pkpc of the line of sight to the QSO. The 
HIRES data are in black, while the red shows our Voigt profile decompo- 
sition of the H I absorption near the redshift of the galaxy. The continuum 
and zero level of the spectrum are shown in dashed and dotted lines respec- 
tively. The systemic redshift of each galaxy is marked by the vertical dashed 
fine at km s"' . Note that the continuum is depressed in some of the spec- 
tral regions surrounding Q1442-BX333, Q1442-MD50, Q0100-BX210, and 
Q144 2-MD84 by a DLA or sub-DLA near the galaxy redshift as described in 



the random-redshift distribution. The more significant excess 
over random of the 300 km s"' version is primarily due to the 
exclusion of unrelated absorbers at large velocity separation; 
however, we note that at Z)tran> 1200 pkpc, the differential 
Hubble velocity associated with this distance along the line of 
sight is > 300 km s"', meaning that for large D^-^n it may be 
more appropriate to adopt max(A^Hb700km s~') as the relevant 
statistic. Regarding the sum statistic for the smaller velocity 
window (not shown), the value of the plateau and that of the 
random sample is ^ 0.3 dex higher for sum(A^Hi,300 km s"') 
than for max(A^Hi,300 km s"'), similar to the variation in the 

statistics shown in Figure |9] 

As we will argue in ^ and ^8.11 the velocity and spatial 
scales of 300 km s"' and 300 kpc capture the most signif- 
icant excess in both the column density and the number of 
absorbers near galaxies. In Figure [TT| we provide the indi- 
vidual measurements max(A^Hi, 300 km s~') for all galaxies 
in the sample with Z)tran< 300 pkpc. Note the large intrinsic 
scatter in max(A^Hi), even at fixed impact parameter For the 
10 galaxies with the smallest impact parameters (Z)tran< 100 
pkpc), the relevant portions of the QSO spectra within ±1000 
km s"' of Lya at Zg^i are reproduced in Figure [T2l 

4.2.1. The large-scale distribution of H I 

We now consider the larger-scale distribution H I around 
galaxies. Unfortunately, the sampling of galaxies with Z)tian> 
2 pMpc is comparatively poor in our sample due to the sur- 
vey geometry of most of the KBSS fields (typically 5.5 x 7.5 
arcmin on the sky). This scale is imposed by the footprint of 
LRIS; with each field roughly centered on the bright QSO, 
the maximum observed impact parameter would be ~ 2.15 
pMpc at (z) =2.3. However, three out of 15 survey fields (see 
Table [T]i were imaged with other instruments covering larger 
angular sizes and thus provide information on larger trans- 
verse scales. In Figures l9l and [TOl we use wider bin sizes for 
absorbers with Z)tran> 1 pMpc in order to consider the large- 
scale distribution. This reduces the shot noise in the bins with 
Dtmn> 2 pMpc. Figures |9] and [To] demonstrate that max(A^Hi) 
remains higher than the global median value in the IGM (dark 
horizontal line) out to ^2 pMpc, and then begins to decline. 
For larger Dnan, the data suggest column densities at or below 
that of random places in the IGM. 

Again, considering the degree of scatter in max(A^Hi) at 
fixed impact parameter, Figure[T3]shows max(A^Hi,700km s"') 
for various bins in Dtran, as indicated. The top row of panels 
correspond to Z)tran< 300 pkpc, while the bottom two rows of 
panels consider larger impact parameters. Each panel shows 
the Kolmogorov-Smirnov probability that the two histograms 
are drawn from the same parent distribution. Notably, only 
the 2 <£>tran< 3 pMpc bins have a distribution of max(A^Hi) 
consistent with the random sample. 

Thus, we have shown that the column density of H I peaks 
sharply at the position of galaxies in the transverse direction, 
that the width of the peak is ~ 300 pkpc, and that there re- 
mains a significant excess of H I gas to Z)tran— 2 pMpc. In ^ 
we discuss the implications of these results. 

4.3. 3D Distribution ofNm 

The 3D distance, D^n, is computed using the quadrature 
sum of the physical impact parameter (Dtran) and the line-of- 
sight distance calculated assuming the velocity differences Av 
are due entirely to the Hubble flow. 
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Figure 13. The maxCWm, 700km s"' ) statistic as a function of impact parameter (different panels). The hatched histograms are the values for the random sample, 
the solid histograms are the values for the real sample. These histograms quantify the variation from galaxy to galaxy of the max(Nii[) statistic at fixed impact 
parameter. Pks is the probability that the two max(A'Hi) sets were drawn from the same distribution. Notably, the last panel with the highest value of Aran has 
the least significant departure from the random distribution. 

decline of A^hi as a function of D3D is quite smooth, but again 
strongly peaked at the position of galaxies. The pixel anal- 
ysis of the KBSS sample recently completed by iRakic et alj 
I^Ollb) studies in detail the 3D distribution of H I optical 
depths. There is excelle nt agr e ement between the optical 
depths measured in Rak ic et alJ (1201 Ibi) and the A^hi trends 
shown in Figure [14] The smoothness of the decline is caused 
by the shifting of absorbers with small Dtran and modest | Av| 
into bins at larger D3D. Thus, the signal appearing in the in- 
ner 300 pkpc in Figure|9]is distributed across a larger number 
of bins in Did as a result of the velocity distribution of the 
absorbers. 

The nature of the velocity width of the excess absorption 
will be discussed at length in ^ however, it should be noted 
that all measurements of |Av| rely on the accuracy of the 
galaxy redshifts and in addition are affected by whatever pe- 
culiar velocities are present, whether due to random motion, 
inflows, or outflows. 



D3D(Av,Z,Aran) ^ 



'(Aran)2 + 



Av 
H(z) 



where H(z) is given by 

H(z) = i/oVf^m(l+z)3 + OA 
such that H(z = 2.3) in our cosmology is 

H(z = 2.3) = 240 km s~' Mpc"' . 



(7) 



(8) 



(9) 



(10) 



In this formalism, each absorber has a unique D^d with re- 
spect to a galaxy in the same field. 

The 3D distance, due to its strong dependence on | Av|, re- 
quires that absorbers have both small Dn-an and very small val- 
ues of I Av'l in order to populate bins at small values of Dt,u- 
As a result, D^u has the effect of isolating those absorbers 
most likely to be associated with the galaxy without imposing 
a velocity cut. 

Figure[T4]shows that the median A^hi stays above that of an 
average place in the universe out to D^u — 3 pMpc.^^ The 



4.4. Connection to Galaxy-Galaxy Pair Results 



Steidel et al. (2010) presented a sample of 512 close angu- 
lar pairs of galaxies with different redshifts (drawn primarily 
from the same KBSS catalogs used in the present paper), us- 
ing the spectrum of the background galaxy to probe gas as- 
sociated with that in the foreground. They were able to mea- 



Note, this is consistent with the quadratic sum of the extent of the ve- 



locity and transverse distance excesses seen in previous sections. 
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Figure 14. The median column density of all absorption systems witliin 
±1400 km s"' of a galaxy as a function of the 3D distance between the 
absorber and the galaxy. The symbols have the same meaning as those in 
Figure|9] Note the steeply rising column densities at small D^jy and that the 
median value remains above that of a random location (horizontal bar) out to 
~3 pMpc. 

sure the strength of absorption from H I and several metal- 
lic species over a range in impact parameter Aran= 8- 125 
pkpc. The principal advantage of this method is that it al- 
lows for probes at very small angular separation {9 < 5")- 
obtaining statistical results for galaxies at such small sepa- 
rations from QSOs is difficult due both to the "glare" of the 
QSO and the relative rarity of QSO-galaxy pairs with very 
small separations (the smallest QSO sightline-galaxy separa- 
tion is Z)tran= 50 pkpc or w 6") . As such, the galaxy pair 
results are highly complementary to the QSO sightline study 
described in this paper 

Using background galaxies instead of QSOs results in two 
important differences between these studies: first, background 
galaxies have a projected "beam size" of 1 kpc at the lo- 
cation of the foreground galaxy, whereas background QSOs 
have a projected beam of order ^ 1 pc. As a result, the ab- 
sorption seen against background galaxies measures a combi- 
nation of the covering fraction of gas on kpc scales and the 
column density of absorbers. Further, background galaxies 
are faint and theref ore only low-disper sion, low-S/N spec- 
tra can be obtained. ISteidel et al.l ( 12010!) thus used stacks of 
background galaxy spectra shifted into the rest frame of the 
foreground galaxies to quantify the average absorption profile 
surrounding these galaxies. The lower-resolution spectra do 
not separate into individual components as would be found in 
individual QSO spectra. As such, the galaxy pair method al- 
lows for the measurement of eq uivalent widths (Wb) only. As 
discussed bv lSteidel et al.l (12010 ) interpretation of Wo is com- 
plicated by its sensitivity to the covering fraction and velocity 
extent of the absorbing material, and its relative insensitivity 
to ionic column density. 

Steidel etjp (12010 ) measured Wo of H I, C IV, C II, Si IV, 
and Si II in bins of impact parameter (Dtran)- The large Wo, 
particularly at Dnan ^ 40 pkpc, seemed to require high line- 
of-sight velocity spread in the gas- likely higher than could 
be easily explained by gravitationally-induced motions, but 



easily accounted for by the high velocities observed "down 
the barrel" in the spectra of LBGs, which typically reach 
Av ~ 800 km s~' for galaxies with properties similar to those 

i n our sample. 

ISteidel et alJ (|2010) found that a simple geometric and kine- 
matic model of outflows from LBGs could account simultane- 
ously for the behavior of Wo as a function of Dtran as well as 
the shape of the blue-shifted line profiles of strong interstel- 
lar absorption features observed along direct sightlines to the 
same galaxies. The presence of detected C II absorption to 
^^tran ^ 90 pkpc was citcd as evidence for a non-negligible 
covering fraction of H I gas h aving A^ht> lO'^ cm ~^ at such 
large galactocentric distances. ISteidel et al.l (120101) remarked 
that 90 kpc is very close to the expected virial radius r^\^ for 
LBGs similar to those in the current sample. Our max(A^Hi) 
statistic agrees well with this inference (Figure|9l); in fact, we 
find a covering fraction ~ 30% for absorbers with A^hi> 10'^^ 
cm"^ for r ^ r^,\r- see §15.21 



In addition, ISteidel et all (1201 Oh used the same HIRES 
spectra as this work to show that the Wo(LyQ;) measured from 
the relevant portions of the high resolution spectra, averaged 
in the same way as the background galaxy spectra, are consis- 
tent with an extrapolation to Dtran — 250 pkpc of the trend seen 
in the gala xy-galaxy pairs resu lts; a similar conclusion was 
reached by iRakic et al] (1201 Ibl) . also using equivalent width 
analysis of the QSO spectra, where a smooth trend was noted 

out to Z)tian> 1 pMpC. 

4.5. Comparison to previous studies at z> 2 

lAdelberger et al.l (120031 12005 a') conducted the first system- 
atic studies of high-z galaxies and their surrounding IGM us- 
ing sightline surveys of the 2 < z < 4 IGM paired with large 
LBG surveys des i gned t o probe galaxies in the same volume. 
lAdelberger et alj (l2003h analyzed the transmitted flux in the 
Lya forest of background QSOs, evaluated near the redshifts 
of survey galaxies. At the time, this was most easily accom- 
plished using z 3 LBGs and z — 3.5 background QSOs. 
These authors did not attempt Voigt profile decompositions 
(as in the present work) but focused on transmitted flux be- 
cause it made for easier comparisons to theory and because 
the spectra covered the Lya transition only, making measure- 
ments of column density or optical depth difficult due to lim- 
ited dynamic range for any A^hi ~ l O'^^^cm"^. Based on 8 
QSO sightlines covering 3 < z ^ 3.6. lAdelbergeret alJ (120031) 
found that excess H I absorption (i.e., lower transmitted flux 
than the average IGM at the same redshift) was present within 
^ 5h"'comoving Mpc (cMpc) of galaxies [1.7 physical Mpc 
(pMpc) at (z) = 3.3 using the cosmology adopted in the 
present work], and an intriguing but not highly significant 
lack of H I very close to these galaxies (< 0.5h"'cMpc or 
< 170 pkpc). C IV absorption was observed to be correlated 
with galaxy positions to out to 2.4h"'cMpc (or 800 physical 
kpc). The cross-correlation of C IV systems with LBGs was 
found to be similar to the LBG auto-correlation, suggesting 
that metal enriched IGM and galaxies shared the same vol- 
umes of space. The strongest C IV absorbers were so strongly 
correlated with LBG positions that the authors concluded that 
they must be causally connected to one another 

Adelbergeret al. (2005a) extended similar studies to 1.8 < 
z < 3.3, using a larger number of QSO sightlines and a wider 
range of QSO spectra for the analysis. Based once again on 
the transmitted flux statistics, the large-scale excess A'hi near 
galaxies was consistent with the results of .Adelberger et all 
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Figure 15. The covering fraction, /c, of absorbers for various Nyii thresliolds (different panels) as a function of -Dtran- The solid histogram represents the fraction 
of galaxies with an absorber of a given Nm or greater, within ±300 km s"' . The hatched histogram represents /c at random locations in the IGM; the horizontal 
dashed line marks /c = 1 (100% covering). The bin size is 100 pkpc for absorbers with Aran < 1 pMpc and 200 pkpc for those with Aran > 1 pMpc. 

(|2P03), but the "turnover" of H I absorption on the smallest 
scales was not confirmed using the larger sample at some- 
what lower redshift. Generally, the transmitted flux was de- 
creasing with galactocentric distance, though it was still the 
case that ^ 30% of the galaxies with the smallest impact pa- 
rameters showed a lack of strong absorption, interpreted as 
evidence that the gas is clumpy. The results for C IV were 
extended, and it was shown that the correlation with galaxies 
grows increasingly strong as A^civ increases; the correlation 
length of A^civ ^ 10'^'^ cm"^ absorption systems was similar 
to that of the autocorrelation length of the galaxies. The au- 
thors attempted to correlate the observed IGM properties with 
galaxy properties, but no significant correlations were found 
given the relatively small sample of galaxies at small impact 
parameters. 

There is no overlap in the data sample used in this paper 
with that of Adelberger et al. (2003) which focused on higher- 
redshift galaxies and QSOs. Three of the fields (Q1623, 
HS1700, and 02343) included in our analysis were also in- 
cluded in lAdelberger et al.l ( l2005ah which considered similar 
redshifts to this work; however, we have increased the S/N of 
the QSO spectra and also added many galaxy redshifts to our 
catalogs for these fiel ds since the earlie r analysis. The most 
surprising result of the lAdelberger et aH (|2003) sample, espe- 
cially in Ught of the work presented here, was the reported 
deficit of absorption found within 0.5 cMpc or 170 pkpc 
of galaxies. Adelberger et al. (2005a), with a larger data set, 
found that 1/3 of galaxies had relatively-little H I absorption 
(consistent with th e pixel statistics of our sample presented in 
iRakic et ani201 Ibh . but that the majority of galaxies were as- 
sociated with significant absorption. In our sample, only 1/21 
galaxies with Z)tran< 170 pkpc has a sum(A^Hi, 300 km s"') 
value less than the median of the random sample. As such we 



concur with the argument of Adelberg er et al.l (j2005a) sug- 
gesting that the reported deficit was due to the small sample 
size presented. Further, we note that a measure of A^hi shows 
that the majority of galaxies do have excess H I absorption in 
their surroundings. 

On larger physica l scales, lAdelberger et al.l fe003) and 
ICrighton et al.l (1201 ll) conside red the large-scal e distribution 
of H I absorbers at z « 3 while lAdelberger et al.l (|2005a) stud- 
ied that at z « 2.5. All found evidence for increased absorp- 
tion to D^D ~ 5-6 cMpc or ~ 2 pMpc. This scale is 
roughly consistent with our measurements of the 3D distri- 
bution of absorption presented in Figure [141 as well as with 
the tren ds in optical depth vs. D^d presented in iRakic et al.l 

(1201 m . 

5. THE COVERING FRACTION AND INCIDENCE OF H I 



As discussed in 34.11 the scale of the strongest correlation 
between A^hi and the positions of galaxies is found within 300 
km s~' of Zgai. Adopting this value as the characteristic veloc- 
ity scale for circumgalactic gas, we can examine other useful 
measures of the gas distribution around galaxies. The cov- 
ering fraction (/c) and the incidence of absorbers (?7abs) are 
two ways of quantifying the geometry of the distribution as a 
function of both impact parameter and A^hi- 

First, we define the covering fraction, /c(£'tran. No), as the 
fraction of galaxies in a bin of impact parameter, Daan, that 
have an absorber within |Av|< 300 km s"' with A^hi > Nq. 
This is equivalent to the geometric fraction of the area of an 
annulus centered on the galaxy that is covered by gas with 
A^Hi>No and |Av|< 300 km s"'. This quantity measures the 
variation within the sample of the decline of A^hi as a function 

of Aran- 

A related quantity is the incidence of absorbers, ryabs, de- 
fined to be the number of absorbers per galaxy within a given 
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Figure 16. The incidence of H I absorbers, rjabs, as a function of impact parameter, Dtran- The solid histogram represents the mean number of absorbers per 
galaxy within ±300 1cm s"' at the given distance, whereas the hatched histogram represents the average incidence of absorbers near randomly-chosen redshifts. 
The different panels show various ranges of A'hi. The bin size is 100 pkpc for absorbers with Z)tran < 1 pMpc and 200 pkpc for those with Duan > 1 pMpc. The 
incidence of absorbers exceeds the random distribution for A'hi > lO'^^ cm"^ to Dtran > 2 pMpc. Note that only absorbers with log(A'Hi) > 14.5 show strong 
association with the positions of galaxies. Table|2]gives the data values for fc and Pe determined from these distributions. 



Table 2 

Absorber Incidence (rjabs) and Excess Probability (Pe) for | Ai'|< 300 km s"'" 



log(/VHi) range'' 




Pe 




Pe 




Pe 


[cm-2] 


Dtnm< 0.1 pMpC 


Aiai,< 0.3 pMpc 


Aran< 2pMpC 


13.0—13.5 


1.1± 0.3 


-0.1±0.3 


1.3±0.2 


O.liO.l 


1.33±0.04 


0.17± 0.04 


13.5—14.0 


1.4± 0.4 


1.0± 0.5 


1.2± 0.2 


0.8± 0.3 


0.87± 0.03 


0.44± 0.06 


14.0 —14.5 


0.6± 0.2 


0.7± 0.7 


0.7±0.1 


1.4± 0.4 


0.53± 0.03 


0.80± 0.09 


14.5 —15.0 


0.8± 0.3 


4.4± 1.9 


0.7±0.1 


3.4± 0.8 


0.32± 0.02 


1.03± 0.13 


15.0—16.0 


0.9± 0.3 


12.7± 4.6 


0.4± 0.1 


4.9± 1.4 


0.22± 0.02 


1.59± 0.20 


16.0—17.0 


0.7± 0.3 


38.2±14.8 


0.2±0.1 


9.7± 3.4 


0.06± 0.01 


1.78± 0.40 


" The values in this table refer to Figure [T6l 



logCWni) range: The range of Whi considered; And the incidence, rjabs. and excess probability. 



Pe considered over three ranges of Dtran- 
assuming Poisson statistics. 



The quoted uncertainties in rjabs and Pe are calculated 



range of A^hi and Aran, and with |Av|< 300 km s"'. Because 
this quantity can be greater than unity, it has larger dynamic 
range and so allows for a more-complete picture of the av- 
erage multiplicity of absorbers at locations close to galaxies. 
Also, because we consider ?7abs in differential bins of A^hi, it 
can be used to measure the degree to which absorbers of a 
given A^Hi associate with galaxies. 

Figure [15] illustrates the dependence of /c on Dttan for var- 
ious thresholds of A^hi (different panels). Within 100 pkpc, 
fc> 0-5 even for A^hi> 10'^ cm~^. Also, every galaxy has an 
absorber with log(A^Hi) > 14.5 within 100 pkpc and |Av|< 
300 km s"'. The distribution within 2.5 pMpc of lower-A^ni 

To illustrate the difference between /c and »7abs, suppose that all ab- 
sorbers have the same Whi- Then for a sample of two galaxies, one with 3 
absorbers and the other with none, /c=0.5 and r;abs= 1-5. 



absorbers is relatively uniform, with /c> 0-5. However, at 
larger threshold A^hi (bottom three panels) we see that /c is 
high ( ^ 0-5) only within ^2 00 pkpc. We will return to our 
measured values of in ^5. 21 and ^5- 31 where we compare our 
measurements to those made at low-z and also to results from 
numerical simulations- 
Figure [16] shows the average incidence of absorbers, yyabs, 
as a function of distance, Du-^„, and ranges of A^hi (different 
panels)- Notably, in all panels (i-e-, at all columns densities) 
the average 77abs is higher than random for Dtran < 2 pMpc- 
For absorbers with A^hi> lO''*^ cm~^ (three bottom panels), 
there is a clear peak at small values of Dtran- No similar peak 
is present in the distributions of absorbers with A^hi< 10'"*^ 
cm"^ (top panels)- Clearly absorbers with A^hi> 10'"*^ cm"^ 
are more tightly correlated with the positions of galaxies than 
absorbers of lower column densities- 
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Figure 17. The incidence of absorbers, ry^bs. with Whi > lO'^ cm"" as a 
function of impact parameter. The solid histogram represents the distribution 
of these high-Wm absorbers for | Av| < 300 km s"'. The hatched histogram 
represents the average incidence of the same absorbers near randomly-chosen 
redshifts. The vertical dotted line indicates the distance at which the inci- 
dence of the real distribution becomes comparable with that of the random 
distribution, Duan ~ 2 pMpc. The dotted horizontal line marks rjjibs = 1- 

These distributions can also be used to determine the excess 
probability (over that of a random place in the IGM) of inter- 
secting an absorber of a given column density, within a range 
of Aran and within some | Av|. The excess probability, Pe, is 
defined through comparison to the random distribution: 

^abs ^abs. random / 1 l \ 

= ■ (tt) 

^abs .random 

Table |2] summarizes the measured values of ?7abs and Pe for 
absorbers binned in A^hi for the velocity window |Av|< 300 
km s"' and three distance cuts: Z)tran< 100 pkpc, Dami< 300 
pkpc, and Z)tran< 2 pMpc. 

5.1. Absorbers with Nni > 10^'^-^ cm~^ 

We have shown above that absorbers with A^hi > 10'^'^ cm"^ 
appear to trace the positions of galaxies in our sample with 
high fidelity. The incidence of absorbers with A^hi > 10'"* 
cm"^ and |Av| < 300 km s"', as shown in Figure [17] nicely 
encapsulates the "shape" of the CGM. Similar to Figure|9]for 
Dunn < 300 pkpc, one sees rising values of r/abs as the galacto- 
centric distance is reduced. From 300 pkpc < Duan < 2 pMpc, 
?7abs reaches a plateau value with ryabs 0.5 [2 pMpc =1.4 h"' 
pMpc « 4.6 h"' cMpc (at z = 2.3)]. For Aran > 2 pMpc, T^abs 
drops to values consistent with the average IGM. 

Recalling the quantity P(Av,Dtran), defined in ^4.11 as the 
probability, per galaxy, of intersecting an absorber at a given 
Av and within the specified range of Dtran, here we consider 
the velocity distribution of absorbers with A^hi> lO'"^'^ cm"^ in 
bins of Dtran as shown in Figure [TS] 

Fitting a Gaussian to the excess absorbers near galaxies 
compared to random gives a deviation (cr(Ai )) of 187 km s"' 
again suggesting that the velocity scale of the excess is Avw 
±300kms-i. 

Integrating the distribution in the first panel of Figure [18] 





(Av> 




A I' window 




Pe 


[pkpc] 


[kms-'] 


[kms-'] 


[kms-'] 






0-300 


4±33 


164± 28 


-350 to 350 


1.70±0.02 


4.1 


300-600 


-74±50 


253± 49 


-550 to 450 


1.07±0.01 


1.19 


600-900 


29±49 


389± 49 


-750 to 850 


1.31±0.01 


0.80 


900-1200 


-43±78 


407± 87 


-850 to 750 


1.60±0.01 


1.14 


1200-1500 


61 ±78 


335± 80 


-650 to 750 


1.20±0.01 


0.76 


1500-1800 


42±96 


344±109 


-650 to 750 


1.06±0.01 


0.63 



" The values in this table refer to Figure [TS] 

Z5tran:The D(ran window considered; (Av): the Gaussian velocity cen- 
troid and its associated error; cr^^,,^ : the standard deviation of the Gaus- 
sian fit and its associated eiTor; Ai' window: the velocity window used to 
compute r^abs and Pe (comparable to ±2(7); rj-^bs- the incidence, Pe: the 
excess probability. The quoted uncertainties in r^abs ^re calculated assum- 
ing Poisson statistics. The error in Pe is dominated by the variation in 
individual realizations of the random sample: typically the variation in Pe 
is ~ 0.05 except for the first bin (0 < Aran< 300 pkpc) where the variation 
is ~ 0.2. 

corresponding to Duan < 300 kpc results in an incidence aver- 
aged over ±350 km s"' (~ 2a) and 300 pkpc, ryabs = 1 -6. Sim- 
ilarly, comparing the value of Tyabs measured in the real and 
random distribution one infers that strong absorbers (Nm> 
lO''*'^) are > 4 times (Pe = 4.1) more likely to be found within 
300 kpc and ±350 km s"' of a galaxy in our sample, than at a 
random place. The parameters of the Gaussian fits in Figure 
[Tsl as well as the inferred r/abs and Pe, are listed in TableQ] 

In this section, we have shown that the covering fraction, 
/c, of absorbers with log(A^Hi) < 14.5 is roughly uniform and 

> 0.5 out to Dtran ~ 2 pMpc, but that those with higher A^hi 
only have /c > 0.5 within Duan < 200 kpc (Figure [TSll. Con- 
sideration of the incidence, T^abs, indicates that absorbers with 
log(A^Hi) > 14.5 are more directly related to galaxies (Figure 
[T6l . Those absorbers nicely encapsulate the "shape" of the 
CGM (Figure [TtT i and are > 4 times more likely to be found 
within |Av| < 350 km s~' and Aran < 300 pkpc of a galaxy 
in our sample than at a random place in the IGM (Figure [TSl). 

5.2. Covering Fractions: Comparison with Simulations 

Motivated by the desire to predict the observational signa- 
tures of cold accretion streams, two recent theoretical papers 
have considered the coveri ng fraction of absorbers of various 
A^Hi surrounding galaxies. iFaucher-Giguere & KeresI ( 1201 Ih 
consider the covering fraction of Lyman Limit (LLS, 17.2 < 
log(A^Hi) < 20.3) and Damped Lyman Alpha (DLA, log(A^Hi) 

> 20.3) absorbers originating within cold streams near two 
simulated star-forming galaxies at z = 2 using cosmological 
zoom-in simulations that do not include galactic winds. Since 
these authors explicitly did not consider galactic winds, their 
covering fractions are approximately^^ a lower limit on the 
expected values. They considered a galaxy slightly less mas- 
sive than those in our sample (Mdm = 10'^^ Mq), as well as 
one comparable to those in this work (Mdm = 10'^ M©). 

For comparison, we consider the fraction of sightlines at 
a given Dtran for which sum(A^Hi,700 km s"') falls within the 

^' These values are only a lower-limit in the case that outflowing winds 
do not affect the gas distribution within the filamentary cold streams, 
van de Voort et al. (2011b) found that the wind prescription in simulations 
has little effect on the amount of gas that accretes onto halos but does 
affect the amount which is delivered into the ISM of the galaxy itself. 
IFaucher-Giguere et alj (201 11) argue that energetic winds can undergo hydro- 
dynamical interactions with cold-streams removing some of the inflowing 
gas. 
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Figure 18. The velocity of absorbers with A'hi> lO'*^ cm ^ in bins of Dtran, as indicated. Note that Nm > 10''" cm ^ absorbers are > 4 times more likely to 
be found within ±300 km s"' and 300 kpc of a galaxy than at a random place in the IGM. See Table[3]for r?abs and the excess probability for all of the panels. 



Table 4 

Fc : Comparison wi th IFaucher-Giguere & KeTeg <20Tll) 



Sample 


logC^Hi) [cm-2] 


Fc«l rvi,-) 


Fc«2 rv,r) 


Mdm = 10'2 M© 


17.2-20.3 
> 20.3 


11% 
4% 




Mdm = 10" 5 Mq 


17.2-20.3 
> 20.3 


12% 
3% 


4% 
1% 


This work" 


17.2-20.3 
> 20.3 


30±I4% 
0+10%'' 


24±9% 
4±4% 



" The Fc tabulated here from "this work" are the fraction of galaxies 
with a sum(A'Hi) statistic in the column density range. This is subtly 
different from the /c presented earlier. Note that while we use a 
velocity cut of |Av|< 700 km s"' here, the results would be the 
same if we used a ±300 km s"' window. Had we considered fc as 
opposed to Fc, we would have calculated 20% for the LLS covering 
within both 1 and 2 Tyi, . The uncertainties quoted for "this work" are 
l-cr en'ors calculated assuming Fc follows a binomial distribution. 
'' We calculate the 1-a upper hmit on this non-detection as 
1 -(1 —0.68)'/''+"' where n = 10 is the number of systems consid- 
ered. 

specifie d column density range. This is most akin to the re- 
sults of Faucher-Giguere & Kerel (1201 li) as they measure the 
column density of absorbers after projecting their simulated 
cube onto a 2D plane. Note that here we use differential bins 
in A/hi(17.2 < log(A/Hi) < 20.3), and since we are considering 
sum(A/Hi) rather than max(Nm), we use to denote the asso- 
ciated covering fraction (to differentiate it from the /c mea- 
sured earlier in this section). 



Our measurements are compared with the 
iFaucher-Giguere & KeresI (1201 ll) simulation results in 
Table ID Particularly in the case of LLS gas, we find ^ 3 
times higher within rvir and ^ 6 times higher within 2rvir 
compared with the simulations. However, if the simulations 
results are treated as lower limits, then clearly they are 
consistent with the observations. 

iFumagalli et al. (2011) also considered the covering frac- 
tion of H I surrounding 6 LBG-type galaxies at z = 2 - 3 us- 
ing cosmological zoom-in simulations with galactic winds in- 
cluded (though they are relatively weak with their particular 
implementation). They consider absorbers within and 2rvir 
for various thresholds in log(A/Hi) : > 15.5, > 17.2, > 19.0, 
and > 20.3 cm~^. A compariso n with the observations is given 
in Table |5] The models of .Fumagalli et al.l (1201 ll) clearly 
under-predict the presence of 15 .5 <log(A^Hi) < 17.2 gas 
surrounding galaxies. Similar to i Faucher-Giguere & KeresI 
(|2011), they also seem to under predict the radial extent of 
gas with A^Hi ^10'^ cm"^. 

5.3. Evolution of the CGM from z 2.3 f o z ^ 1 

There have been numerous studies of the low-z galaxy-lGM 
connection. In the interest of brevity, below we compare our 
measured fc with a small number of studies with comparable 
size and statistical power. 

As discussed bv lDave et al] ( 119991) and lSchavd (•2001), due 
to the expansion of the universe, the collapse of structure, 
and the evolution in the intensity of the metagalactic ioniz- 
ing background, the A^hi associated with a fixed overdensity 
declines with redshift (e.g., an absorber at z = would have 
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Figure 19. The covering fraction of high- and low-redshift absorbers in cumulative bins of Diran (i.e., at each value of Aran, Do, we consider all absorbers with 
smaller Aran, as opposed to those absorbers with Do < Aran < Di.) In the (light-blue) solid histogram we plot our high-z data for /c close to galaxies. The 
(dark-blue) cross-hatched histogram shows the distribution at random places in the high-z IGM. The (red) vertically-hatched histogram shows the distribution of 
fc of low-z absorbers taken from lWakker & Savagel 42002). Here we consider the same cut in Nm at both epochs (Wo = 50mA is equivalent to Nm = lO'^ cm"^ 
and similarly Wo = 300mA is equivalent to A'hi = lO''' cm"^). 



Table 5 

Covering Fraction: Comparison with lFumagalli et aH 1201 ID 



Sample 



log(A'Hi) [cm 



fc«l/-vir) /t«2r„r) 



Fumagalli et al. 


> 15.5 

> 17.2 

> 19.0 

> 20.3 


38% 
16% 
6% 
3% 


22% 
7% 
3% 
1% 


This work" 


> 15.5 

> 17.2 

> 19.0 
>20.3 


90±9% 
30±14% 
10±9% 
0+'o% 


68±9% 
28±9% 
8±5% 
4±4% 



" Values in "this work" are computed as described in Table|4] Were 
we to consider | Ai'| < 300 km s"' rather than | Ai'| < 700 km s"' , 
only the Fc for log(A'Hi) > 15.5 within 2 ry\i would change. The 
value for 300 km s"' would be 60±10%. 

A^Hi ^ 20 times lower than would be measured at the same 
overdensity at z ss 2.3). As a consequence, absorbers with 
the same A^hi at different redshifts trace different structures, 
and have very different incidence rates. For example, there 
are no regions of the Lya forest at z ^ 2.3 that have zero 
H I absorbers in a velocity window of ±700 km s"' - ev- 
ery galaxy in our sample can be associated with (generally) 
a large number of absorbers. This is not the case at low-z, 
where absorption-line spectra are sparsely populated and gen- 
erally galaxies would be associated with very few (possibly 
no) H I absorbers. Similarly, the continued growth of struc- 
ture and the rapid decline in the star-formation rate density 
toward z ^ means that a typical low-z galaxy is expected 
to be in a very different place on its evolutionary sequence 
compared to a "typical" galaxy at z > 2. Large differences are 
also expected for gas fractions and baryonic accretion rates at 
low-z compared to high. Thus, one might reasonably expect 
circumgalactic gas to differ as well. 

With these caveats in mind, we compare the covering frac- 
tion, /c, of the (z) = 2.3 sample with results at low-z. 

5.3.1. The CGM at OA < z < 0.9 



IChen et alj (1200 Ih considered 47 galaxies with (z) = 0.36 
and Z)tran< 330 pkpc Their galaxy sample covered a wide 
range of properties, with 68% having Lb > 0.25 Lg. The QSO 
spectra used were generally not of sufficient quality to allow 
for Voigt profile decompositions and thus the au thors mea- 
sured t he equivalent width (Wo) of Lya absorption. iChen et alJ 
(120011) considered a bsorbers to be associated with galaxies if 
I Av| < 500 km s"' (iMorris & Jannuzi 2006). 

For Wo > 350 mA (A^hi ^ 10'"^ cm'^) and Dtra„< 330 pkpc, 
34/47 galaxies in their sample had detectable Lya absorption, 
corresponding to fc= 0.72. In the (z) = 2.3 sample there are 
48 galaxies with Z)tran< 330 pkpc. Using the same velocity 
window |Av|< 500 km s"', we find 40/48 (83%) have Lya 
absorption withA^HT> lO'^ cm ~^, consistent with the measure- 
ments from Chen et alJ (1200 Ih . For clarity, the results from the 
low-z studies as well as our measurements are reproduced in 
Tabled 

5.3.2. The CGM al 0.005 < z < 0.4 

iProchaska et alj (1201 Ih studied the association of z < 0.2 
galaxies i n 14 fields s urrou nding background QSOs. Simi- 
lar to the IChen et aTl (120011) study, they generally relied on 
conversions of published Wo measurements into A^hi using an 
assumed bd- In total, their galaxy survey included 37 galaxies 
with Aran< 300 pkpc and 1200 galaxies in all with (z) = 0. 18. 
They use a velocity window of |Av|< 400 km s"'. For the 26 
galaxies with L > 0. 1 L* within 300 pkpc, 25 had accompa- 
nying Lya absorption with A^hi> lO'"* cm"^, or fc=96%. In 
our sample, there are 43 galaxies with Z)tran< 300 pkpc, all 
of which have a A^hi> lO'"' cm"^ absorber within |Av|< 400 
kms"' (/c=100%). 

These values appear consistent; however, the covering frac- 
tion of A^Hi> lO'^* cm"^ gas is a very poor measure of the ex- 
tent of the CGM at high redshift. For example, if we consider 

The Dtra n rang es quoted here have been adjusted from the cosmology 
considered in IChen e t al. ( 2001) to the cosmology used in this paper at the 
mean redshift of the sample, (z) = 0.36. 
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Figure 20. Same as Figure [19] but here we consider different A'hi cuts for the high-z sample. We have found the A^hi lower limit for the high-z sample which 
best reproduces the trend of the low-z /c as a function of Dtran ■ We can see that the trends of high-z absorbers of 1 .5 dex higher Nm are well matched to the 
distributions of low-z absorbers. 



the same velocity window 



A^Hi threshold (A^hi> 10'^ cm-^ 



') and the same 
= 95% within 1 



Av|<400kms- 
), we find /c 

pMpc and /c= 94% within 2 pMpc; i.e., there is no apprecia- 
ble change in the fc of g as of this A^hi with increasing distance. 
IProchaska et al. dfoi lb find /c=70% for Dtran < 1 pMpc in 
their sample. 

Measurements with a higher A^hi threshold have larger dy- 
namic range ; for A^Hi> lO'"* cm'^, IProchaska et"an (l2011h 
measure fc=lQ% for Aran< 300 pkpc, and /c=38% for 1 
pMpc. For our (z) ^ 2.3 sample, we measure f^= 81% for 
Aran <300 pkpc and /c=70% for Aran <1 pMpc. 

5.3.3. The CGM in the Local Universe: cz < 6000 km 



IWakker & Savaj^ (l2009l) (hereafter, WS09) studied the 
gaseous distribution around local galaxies {cz < 6000 km s"') 
along 76 lines of sight. The galaxy catalog considered in- 
cluded ~ 20,000 local galaxies. The absorber catalog in- 
cluded 115 intergalactic Lya absorbers (i.e., not from the 
Galaxy) over the same redshift range. They calculate /c in 
cumulative distance bins (e.g. Dtran < 200 pkpc). 

For comparison to the WS09 sample, we use a velocity win- 
dow |Av|< 400 km s"' with respect to galaxy redshifts in 
our high redshift sample. From WS09 (Table 10) we use the 
values of fc measured for the galaxy sample with L > 0.25 
L*. They consider two classes of absorbers, those with 
Wo > 50mA (A^Hi> 10'^ cm'^) and Wo > 300 mA (A^hi ^ 10'"* 
cm"^). We make a direct comparison to our sample with the 
same Nm threshold in Figure [19] and Table |6] Our measure- 
ments are shown in the solid histogram, the random high-z 
sample is plotted in the (dark-blue) cross-hatched histogram, 
and the WS09 sample is represented by the (red) vertically- 
hatched histogram. One can clearly see that at all Dtran and for 
both thresholds in A^hi, the z ^ sample shows significantly 
less covering than the high-z sample. 

As mentioned above, absorbers of fixed A^hi are theoreti- 
cally expected to trace higher-overdensity gas at low redshift 
than at high redshift. By considering various A^hi thresh- 



Our spectroscopic sample includes galaxies to M - 
alent to L > 0.25 L* at z 2.3. 



25.5. This is equiv- 



olds in our own data, we found that /c of the low-z W) > 
50mA(A^Hi> lO'"' cm~^) sample was most comparable to high- 
z absorbers with A^hi> lO'"*'^ cm"^; similarly, the WS09 Wo > 
300mA(A^Hi ^ 10'"* cm"^) can be matched to our high-z ab- 
sorbers with A^Hi> 10'^'^ cm~^. Note that we chose the 
A^Hi threshold which best reproduced /c measured on ~^Mpc 
scales. 

These two comparisons, shown in Figure |20l exhibit re- 
markably similar patterns in the covering fraction of A^hi as 
a function of Dtran- It therefore appears that low-z absorbers 
with log(A^Hi) = No on average trace the same physical re- 
gions around galaxies as gas with log(A^Hi) = No + 1 .5 at 
(z) = 2.3. While we do not directly measure universal over- 
density, these results lend observational support to the pre- 
dictions by Dave et al. ( 1999) and Schaye (2001) that gas at 
the same universal overdensity at low- and high-z would have 
lower A^Hi at low-z. These measurements also suggest that the 
physical size of the CGM around a typical galaxy is 300 pkpc 
at both z and (z) =2.3. 



6. MAPPING THE CIRCUMGALACTIC MEDIUM 

Another way to visualize the distribution of neutral hydro- 
gen surrounding star-forming galaxies is to "map" the distri- 
bution along the transverse direction and the line of sight at 
the same time. As in ^4.31 we use |Av| to compute a line-of- 
sight distance assuming it is entirely due to the Hubble flow. 
In Figures |2T] - |24l the vertical axes represent the Deubbie ve- 
locity scale line-of-sight distribution of gas. The horizontal 
axes correspond to the physical impact parameter between the 
galaxy and the QSO line of sight, Dtran- Colors encode the in- 
cidence, 77ab,s, of the gas in a given pixel in the map and the 
black and white shading represents the signal-to-noise ratio 
(S/N) with which 7yab,s is detected. 

The S/N is determined using Poisson statistics. For a given 
bin in Dtran and Deubbie, we count the total number of absorp- 
tion lines found in the random distribution, Wabs.ran, and the 
number of random locations considered at that impact param- 
eter, Hgai.ran- We also couut the number of real galaxies at that 
impact parameter, Wgai- The "noise" level of the map is then 
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Figure 21. The incidence of 14.5 < log(A'Hi) < IV absorbers. Far Left: Map of the 7?abs with respect to the positions of galaxies (at the origin). Second from 
Right: The distribution at random locations in the IGM. The greyscale panels show the S/N ratios of the maps calculated as discussed in g6] The "bin" size 
for this map is 200 x 200 pkpc and the smoothing scale is 200 kpc in Duan and 140 km s"' (600 pkpc) along the line of sight. The FWHM of the gaussian 
smoothing beam is represented by the elHpse in the top right-hand comer of the 3rd panel from the left. Note that there is a strong peak in the incidence of high 
column density absorbers at close galacto-centric distances. The color bars in the ijabs maps have been re-normalized by Xv = 6-3 in order to match the values of 
T^abs summed over 1 Ai'l < 300 km s"' . The vertical structure in the S/N maps is caused by the variation in the number of galaxies in each bin of Arm- 
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Figure 22. Same as Figure|2T]but zoomed in to show the small scale distribution. The "bin" size for this map is 100 x 100 pkpc. The smoothing scale is 100 
kpc in impact parameter and 140 km s"' (600 pkpc) along the Une of sight. The color bars for the incidence maps have been re-normalized Xv = 12 6 in order to 
match the »7abs summed over [ Av|< 300 km s"' . 



taken to be: 

A/^gal (^abs,ian/^gal,ian) 

<^sis = ^ , (12) 

«gal 

where the quantity inside the square root is the number of ab- 
sorption systems expected (based on the incidence for the ran- 
dom sample) given the number of galaxies in the real sample 
at a given Dii-^„. The square root of this quantity represents the 
Poisson uncertainty in the number of absorbers that would be 
detected; the division by «gai results in an expression for the 
error in r/abs per galaxy. This value is akin to the shaded error 
bars surrounding the median value of the random sample in 



Figure|9] The signal-to-noise ratio is then taken to be: 

S/N = — (13) 

Csig 

The maps are generated using three different bin sizes. We 
consider the intermediate bin size (200 x 200 pkpc) for two 
cuts in A^Hi (Figures 1211 and |23]) which cover impact parame- 
ters out to 3 pMpc. For high-A^Hi absorbers which exhibit a 
peak at small scales, we also consider maps with 100 x 100 
pkpc bins (Figure l22l). For Iow-A^hi absorbers, we provide a 
map with 400 x 400 pkpc bins to emphasize the large scale 
distribution (Figure l24b. After the number of absorbers per 
bin is measured, the values are placed onto a 10 times finer 
grid and smoothed by a gaussian kernel with a full width half 
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Figure 23. Same as Figure llTl but for low column density absorbers with 13 < logC^Hi) < 14.5. Note, that the incidence of these absorption systems is higher 
in the large-scale regions surrounding galaxies, similar to that of high column density systems. However, unlike the high column density absorbers, there is no 
strong peak in the incidence rate of Iow-A'hi systems at small galacto-centric distances. The color bars for the incidence maps have been re-normalized Xv = 6.3 
in order to match the jj^bs summed over | Ai'| < 300 km s"' . 
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Figure 24. Same as Figure l23l but with 400 x 400 pkpc binning. Note that the line-of-sight distribution is compressed compared to the transverse distribution. 
The color bars for the incidence maps have been re-normalized Xv = 3.15 in order to match the r^abs summed over | Ai'|< 300 km s"' . 



maximum equal to the bin size in impact parameter (x-axis) 
and a velocity scale of 140 km s"' or 600 pkpc in Hubble 
distance (y-axis). Binning and smoothing is required to make 
maps of this type as the values of impact parameter and veloc- 
ity are discrete. The scale of the smoothing, however, is moti- 
vated by the data. The smoothing in Da-.^„ is selected to obtain 
a large enough sampling of the inner and outer bins such that 
the noise in the random distribution is reduced. The velocity 
scale smoothing correspon ds ro ughly to the amplitude of the 
redshift uncertainties (see ^2.2\ . 

In Figures |2T|-|24] the normalization of the color bar is re- 
scaled from the raw ?7abs per 100 x 100, 200 x 200, or 400 x 
400 pkpc bin by the multipUcative factor Xv, where: 



^Hubble (300 km s ') 1 .26 pMpc 



£>Hubbie (bin) 



£>Hubbie (bin) 



where DHubbie(bin) is either 100, 200, or 400 pkpc. For the 100 
pkpc bins, Xi = 12.6. This renormalization brings the scale of 
77abs into agreem ent with that shown in Figures [l6land[T7l 

Recently, Rakic et al.l (|201 l b) constructed similar maps to 
those presented here using pixel optical depth analysis (in- 
stead of Voigt profile fits) for a subset of the data presented 
here. Many of the conclusions discussed in the follow- 
ing se ction are c orroborated by the independent analysis of 
.RakicetanJIorib. §4). 

6.1. Absorbers with Nni > 10'"^-^ cm~^ 

Shown in Figure|2T|is a map of Tyabs for 10'"*^ < A^hi < 10'^ 
cm"^. Here we compare a map of the absorber distribution 
around galaxies (left-most panel) to a map which represents 
the general IGM as measured using the random distribution 
(third panel from left). Clearly, the incidence of A^hi> 10'"*^ 
(14) cm"^ absorbers is higher near galaxies than in the general 
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Table 6 

H I Covering Fraction: Comparison with Low-^ Studies^ 



Sample 


WoiLya) 
mk 


logC^Hl) 
[cm-2] 


^tran 

pkpc 


|Av| 
km s ' 




Chen et al. 
0.1 <z<0.9 


350 


Z 14 


<330 


500 


72% 


This work 
z~2.3 




>14 


<330 


500 


83±5% 


Prochaska et al. 
0.005 < z < 0.4 


50 
300 
50 
300 


Z 13 
^ 14 
^13 
^14 


<300 
<300 
<1000 
<1000 


400 
400 
400 
400 


96% 
70% 
70% 
38% 


This work 
z~2.3 




>13 
>14 
>13 
>14 


<300 
<300 
<1000 
<1000 


400 
400 
400 
400 


100 3 % 
81 ±6% 
95±1% 
70±3% 


Wakker & Savage 


50 
50 
50 
50 

50 


i 13 
^13 
^13 
^13 
Z 13 


<200 
<400 
<1000 
<2000 
<3000 


400 
400 
400 
400 
400 


75% 
81% 
49% 
48% 
39% 



z< 0.017 



300 
300 
300 
300 
300 



, 14 
, 14 
14 
14 
14 



<200 
<400 
<1000 
<2000 
<3000 



400 
400 
400 
400 
400 



This work 
z~2.3 



46% 
44% 
24% 
22% 
20% 



> 13 


< 200 


400 


100^%"^ 


> 13 


<400 


400 


98±2% 


> 13 


< 1000 


400 


95±1% 


> 13 


< 2000 


400 


94±1% 


> 13 


< 3000 


400 


94±1% 


> 14 


< 200 


400 


92±5% 


> 14 


<400 


400 


82±5% 


> 14 


< 1000 


400 


71±3% 


> 14 


< 2000 


400 


64±2% 


> 14 


< 3000 


400 


61 ±2% 



" The distance binning is cumulative, so /c is the covering fraction of gas with 
Dinu < the distance listed in the table and |Av| < the velocity listed in the 
table. The low-z samples generally do not include Voigt profile fits, but rather 
measure equivalent widths, Wq. The approximate value of Whi associated with 
each Wo limit is listed. The uncertainties quoted for "this work" are l-cr en'ors 
calculated assuming /c follows a binomial distribution. 

We calculate the I-ct lower limit on this 100% detection as 
(1—0.68)'/"+"' where n = 43 is the number of systems considered. 

We calculate the \-cr lower limit on this 100% detection as 
(1—0.68)'/"+"' where « = 25 is the number of systems considered. 



IGM. 77ab.s remains significantly and consistently higher than 
that seen in the random map within 300 km s"' along the line 
of sight and out to 2 pMpc in impact parameter 

The maps presented in Figure |22] which show the inner dis- 
tance and velocity ranges at higher resolution, allow one to 
compare the scales of the absorption excess in Hubble dis- 
tance and in Dtran in order to detect redshift anisotropics in- 
dicative of particular kinematic patterns. Galaxy redshift er- 
rors and peculiar velocities of gas with respect to galaxies 
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Figure 25. Maps of the median Nm with respect to the positions of galaxies. 
The "bin" size for this map is 200 x 200 kpc. The smoothing scale is 200 
kpc in impact parameter and 140 km s"' (600 pkpc) along the line of sight. 
Left: The map of the true distribution of column densities. Right: The map of 
the column density distribution in the random sample. The significance with 
which we detect structures appearing in this map is quantified with the S/N 
maps presented in the preceding figures. 

will generally expand the line-of-sight velocity distribution. 
Therefore, if the signal along the line of sight is elongated 
more than expected given the redshift errors and the scale 
of the transverse distribution, the remaining velocity structure 
may be attributed to peculiar velocities of the gas with respect 
to galaxies. 

Based on Figure |22] adopting the green contours as rep- 
resentative of the most extreme portion of the overdensity 
of gas (conservatively) corresponds to a velocity scale of 
^300 km s"' and a physical impact parameter ^125 pkpc-". 
Subtracting in quadrature the 140 km s"' smoothing (which 
corresponds roughly to our redshift errors) leaves a residual 
velocity scale of | Av|~ 265 km s"' . This velocity, were it due 
to pure Hubble flow, would indicate a distance of ~' 1 . 1 pMpc. 
In order to consider the portion of the line-of-sight elongation 
likely due to peculiar velocities, we can subtract in quadrature 
the extent of the transverse distribution - in effect subtracting 
a possible Hubble flow broadening. Subtracting the DuubUe = 
125 pkpc or ^ 35 km s"' leaves ~ 260 km s"' ; in other words, 
peculiar velocities of ±260 km s~' are indicated by the data.-*^ 

6.2. Absorbers with Nm < lO'"*-^ cm~^ 

Figure |23] displays a map of the incidence of lower-column 
density absorbers (lO'^*" < A^hi< lO'"*'^ cm"^). Again com- 
paring the left-most panel (galaxy positions) to the third panel 
from the left (random IGM), one sees that the incidence of 
these absorbers is higher in the regions near galaxies than 

One exception to this "rule" is the case of inflowing gas on large scales 
where the inflow velocities counter the relative Hubble flow causing a reduc- 
tion in the hne-of-sight velocity distribution. 

" Maps produced using a symmetric smoothing kernel of 200 pkpc in 
both the hne-of-sight and transverse distance (not shown) also yield excess 
absorption to 300 km s"' along the line of sight. The extent of the line-of- 
sight distribution is robust to reductions in the hne-of-sight smoothing ker- 
nel; however, the S/N associated with the excess is improved by the adopted 
smoothing kernel. 

^- Or, we have seriously underestimated the magnitude of our galaxy red- 
shift errors. 
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in the general IGM; however, the Iow-A^hi map exhibits no 
strong peak at small galactocentric distances. 

Since there is little structure on small scales in Figure l23l 
the map has been more heavily smoothed to emphasize larger 
scales. Drawing attention to the edge of the red contours in 
Figure |24l it appears that the scale of the excess absorption is 
marginally smaller along the DuubUe (line-of-sight) axis (1.5 
pMpc) than along the transverse distance (2 pMpc). This sug- 
gests there may be coherent infall motion on ^ Mpc scales 
compressing the distribution along the line of sight, the Kaiser' 
(1987) effect. Rakic et al. (2011b) analyze the pixel statis- 
tics in the KBSS survey and discuss the kinematics of the cir- 
cumgalactic gas in detail. They find a significant detection of 
compression of the signal along the line of sight which they 
interpret as large-scale infall. 



6.3. Median Column Density Maps 

An alternative visualization of the distribution of H I sur- 
rounding galaxies, avoiding the division of absorbers into low 
and high column density, is a map of the median A^hi in a 
given bin. Figure |25] shows maps where the color bar repre- 
sents the median A^hi in a given bin of Dtran and Z^Hubbie- On 
small scales, the behavior is similar to the high-A^ni absorbers, 
with significant elongation in the line of sight direction, due 
to a combination of peculiar motions of the gas and error in 

Zgal- 

In this section we have shown that the kinematics of ab- 
sorbers with log(A^Hi)>14.5 are consistent with a peculiar ve- 
locity component of ^ ±260 km s"' within ^ 100 pkpc of 
galaxies, and separately that absorbers with log(A^Hi)<14.5 
may exhibit the kinematic signature of retarded Hubble flow 
on scales ^ 1 pMpc. 



6.4. Explaining the Gas-Phase Kinematics 

The kinematics of the gas surrounding galaxies may pro- 
vide important clues about the physical processes occurring 
near galaxies. The two most plausible scenarios for the origin 
of the observed gas is that it traces the large scale structure 
(in which case it is likely falling onto the halo) or that it is the 
result of galactic winds, known to operate in these galaxies. 
Most likely, it is some combination of the two phenomena. 
Unfortunately, the sign of the velocity offset (red-shifted or 
blue-shifted) cannot tell us what the relative motion of the 
gas is with respect to the galaxy since its position along the 
line of sight is not known. For example, gas which is red- 
shifted with respect to the galaxy might be behind the galaxy, 
and thus be either outflowing or moving with the Hubble flow 
at larger distance. Conversely, redshifted absorbers could be 
foreground gas that is redshifted because it is falling onto the 
galaxy. There is no way to know a priori which scenario holds 
in which cases. 

We recall that in ^4.11 it was shown that the full extent in 
Av spanned by absorbers comprising the net excess over ran- 
dom is |Av| ^ 700 km s"', with the most significant portion 
of the excess within |Av'| ^ 300 km s"'. Our map of high- 
A^Hi systems confirms the 300 km s"' excess and allows us to 
measure the extent of the peculiar velocities. After account- 
ing for our estimated galaxy redshift errors (§ 16. Il l, the range 
of significant excess becomes | Av| ^ 260 km s"' . 

In-falling cool gas is likely to be moving at rela- 
tive velocities smaller than the galaxy circular velocity 



dFaucher-Giguere & Kerell201 Ih : 



GMdm , -1 
= 220 km s , 



(15) 



where we assume Mdm = 10'^ Mq and rvir = 91 pkpc (see ^2.11 
and ^8.31 ). Since we measure only the line-of-sight component 
of this velocity, |vios| ^ 200km s"' might serve as an approxi- 
mate upper limit on the line-of-sight velocity component due 
to graviationally induced peculiar velocities. Unless we have 
significantly underestimated the measurement uncertainties in 
Zgai, accreting gas would be expected to have a "quieter" ve- 
locity field than observed. 

On the other hand, outflows observed in the galaxy spectra 
regularly show blue-shifted velocities as high as 800 km s"', 
and large velocity widths appear to be required to explain the 
strength of strongly saturated absorption at modest impact pa- 
rameters (Dti-an ^ 100 pk pc), based on the galaxy-galaxy pairs 
analysis of Stei del et alj (l2010i) . In the context of their simple 
model, the envelope of |Av| is simply the component of the 
outflow speed Voui projected along an observers line of sight; 
this maximum velocity also dictates how strong the resulting 
absorption features will be, with Wo increasing proportion- 
ally to Vout- To model the behavior of Wo with impact param- 
eter, Vout was the principle normalization factor and the in- 
ferred values were 650 ^ Vout ^ 820 km s"' depending on the 
ion. In the present work, the first bin in /Jtran extends from 
50-100 pkpc, with a median D,fan = 8 pkpc, with median 
log(A^Hi) — 16.5. Assuming the lSteidel et alJ (120101) geomet- 
ric/kinematic model (with R^s = 90 pkpc) would predict that 
a line of sight with Z)tian= 80 kpc would have | Av| ^ 0.46vout 
(equivalent to an observable velocity range between ±325 and 
±400 km s"' ) and a threshold A^ei ^ 10'^ cm"^ [analogous to 
our max(A^Hi) statistic], both compatible with the observations 
from the QSO sightlines discussed above. 

7. THE DOPPLER WIDTH 

In addition to A^hi and Zabs, the third component of a Voigt 
profile fit is the Doppler width, b^. In the case of an unsatu- 
rated line, the Voigt profile is well approximated by a Gaus- 
sian in optical depth, with a Doppler parameter: 



FWHM 



(16) 



In the case of purely thermal broadening, the gas temperature 
r can be inferred directly from the Doppler parameter: 



2kT 



(17) 



where m is the mass of the ion and k is the Boltzmann con- 
stant. For H I, 



r=^/,/=4xl04K 
2k 



26 km s 



(18) 



noting that 26 km s"' 
absorber catalog. 

The physical state of 
erned by the balance of 
batic cooling caused by 
photoionization heating. 
In the regime where this 
onic overdensity roughly 



is the median value of bd in the full 



the IGM is predicted to be gov- 
two principle processes: the adia- 
the expansion of the Universe and 
generally from the UV background, 
holds and in the case that the bary- 
traces the dark matter overdensity. 
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Figure 26. The Doppler width, foj, of absorbers with logCWm) > 13 and 
Av|< 700 km s"' as a function of transverse distance from a galaxy. The 
symbols have the same meaning as those in Figure [9] except they refer to 
instead of A'hi- Note the peak in the second bin con'esponding to absorbers 
with 100 < Dcran < 200 pkpc. 

a natural consequence is a relationship between the temper- 
ature of the gas, T, and its density, p. Higher-density re- 
gions, having more gravitational resistance to the Hubble 
flow, cool less and thus have higher temperatures on average 
(see e.g., Hui & Gnedin 1997; Schaye et al. 1999). And in- 
deed, observed distributions of thermally broadened absorp- 
tion lines ex hibit higher for absorbers with higher A^hi 
dPettini et all [19 90; Schaveetal. 2000; Brvan & Machacekl 
120001; lRicottieral.ii2000i: McDonald et al.i i2001). 

Other processes can also broaden individual absorption 
components. For the most diffuse and physically extended ab- 
sorbers, generally with lower A^hi, the Hubble flow itself can 
contribute significantly to the line width. Turbulence (here 
meaning bulk motions of the gas) can also broaden absorp- 
tion features such that: 



2kT 



^turb 



(19) 



where /^turb is the turbulent component of the line width. 

In this section we discuss the distribution of Doppler widths 
(b^) observed in absorption systems as a function of their 
galactocentric distance. The observed trends are much more 
subtle than those in A^hi, but their utility in discerning the 
physical state of the gas motivates a careful analysis. 

7.1. The Dependence ofb^ on Proximity to Galaxies, Dt,an 
Figure |26] shows the Doppler width ibi) of absorbers with 
A^Hi > 10''' cm"^ versus Dt,an- While considerably more 
noisy than the column density trends in the previous section, 
it appears that within the second bin (corresponding to 100 
< £'tran< 200 pkpc) the median value of b^ is considerably 
higher than that of the random sample. 

We consider the trend in b^ versus A^hi for both the sample 
of absorbers close to galaxies and that of the full absorber cat- 
alog. In this way, we can examine the b^ of absorbers at fixed 
Nhi, as shown in Figure |27] in order to evaluate the effect of 
the T - p relationship on this result. Here, each of the four 



closest bins in Dtran (Figurel26]l is broken into individual pan- 
els. The asterisks represent the absorbers close to galaxies. 
The black horizontal lines and shading refer to the full ab- 
sorber catalog. Considering first the top right-hand panel 
with 100 < Z)tran< 200 pkpc (the bin with the majority of the 
signal in Figure l26Tl for A^hi> lO''* cm"-, the median /jj of ab- 
sorbers close to galaxies is larger than the median of the full 
absorber sample at fixed A^hi- This suggests that the larger b^ 
cannot be attributed solely to a dependence on A^hi as in the 
T - p relation. 

Figure [27] shows that the median bi for the sample of ab- 
sorbers within Dtran ^ 100-300 pkpc of galaxies is systemat- 
ically larger than that of absorbers from the full absorber cata- 
logue across most bins in A^hi- However, the absorbers closest 
to galaxies (Aran< 100 pkpc; top left panel of Figure l27Ti do 
not appear to exhibit the same systematic "excess" in median 
bi (cf. Figure |26]|. There are too few galaxies in this inner bin 
to determine whether the result is statistically significant. 

In order to quantify the statistical significance of the appar- 
ent excess b^ in the range Duan < 300 pkpc, we performed 
a Monte Carlo simulation. This was done by drawing ran- 
dom samples of absorbers from the full absorber catalogue 
with the same distribution in A^hi and of the same size as the 
absorber sample close to galaxies. The figure of merit was 
taken to be the average difference computed as follows: in 
each bin of A^hi the difference between the median bj.MC from 
the Monte Carlo sample and the median b^ in the full absorber 
catalog was computed (asterisks minus black bars in Figure 
|27] |. The differences from the A^hi bins were then averaged. 
The fraction of Monte Carlo samples whose average b^ excess 
is greater than that of the Duan < 300 pkpc sightline sample is 
0.005, meaning the excess is significant at approximately the 
3ct level. 

7.2. Doppler Widths vs. 3D Distance 

In Figure|28]we consider the distribution of bd as a function 
of Did, which allows one to consider trends as a function of 
distance without pre-selecting a velocity window. Here again, 
bi is elevated for bins close to galaxies (D3D ^ 400 pkpc). 
In the first bin, the median b^ w 35 km s"' compared to 
bd=27 km s"' for the random sample. 

We again perform a comparison at fixed A^hi and since the 
median bd monotonically decrease with increasing Z)3d, we 
consider cumulative distance bins. For absorbers with D^uK 
600 pkpc, the median bd for absorbers with A^hi< 10'* cm"^ is 
larger than the median measured in the full absorber catalog. 
Notably, for 15 < log(A^Hi)< 16, absorbers close to galaxies 
(D3D< 400 pkpc) have a median b^ > 40 km s"' . If interpreted 
as the result of an increase in gas temperature, the change in 
bd is equivalent to more than doubling T. We also note that 
these intermediate A^hi systems that are characteristic of the 
CGM tend also to be associated with the strongest high ion- 
ization metals (O VI, C IV), often exhibiting evidence that the 
ionization has been produced by shocks in addition to pho- 
toionization (Simcoe et al. 2002). 

As previously stated, the 3D distance is strongly influenced 
by the magnitude of |Av|. Notably, the first three bins in D3D 
(those with the most significant deviation in bd) correspond to 
absorbers with |Av| < 150 km s"'. 

Note that within the full catalog for absorbers whose growth is expected 
to remain approximately linear (A'hi< lO" cm"-, p/p < 10), the median 
values (black horizontal lines) exhibit increasing fej with increasing A'hi, as 
one would expect for the T-p relation. 
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Figure 27. Comparison at fixed Whi of absorbers close to galaxies with absorbers in the full H I catalog. The different panels show four bins of Dam- Asterisks 
represent the median value of in a bin of Nm for absorbers close to galaxies. The horizontal black Unes show the median of absorbers from the full H I 
catalog, and the light shaded boxes represent their dispersion calculated thi'ough bootstrap samples of the same size as those in the real sample. The velocity 
window is | Ai'|< 700 km s"'. In the bins without asterisks, there were no absorbers which satisfied both the Dtmn and N m cr iterion; similarly, asterisks without 
error bars represent bins containing only a single absorber satisflying both the Aran and Nm criterion. Recall from Figure l26l that the strongest signal occurred at 
100 < Aran < 200 pkpc (shown in this Figure in the top right panel). Note that the median fej of all absorbers close to galaxies with A'hi> 10 
than those in the full absorber catalog. 
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Figure 28. Same as Figure l26l but versus the 3D distance. 

In summary, considering the trends illustrated in Figures l26l 
- |29l it appears that absorbers of all A^hi with 100 < Dnan < 
300 pkpc and |Av| < 150 km s"' have significantly larger 
Doppler widths than absorbers of the same A^hi in the full ab- 
sorber sample. 



7.3. Doppler Parameter: Possible Physical Explanations 

The trends in Doppler widths, b^, could have a variety of 
physical origins. As previously discussed, elevated bd could 
be related to the observed trend in A^hi through the T - p re- 
lation in the IGM. However, we showed that at fixed Nui, the 
median b^ of absorbers close to galaxies is still significantly 
larger than in the full absorber catalog (Figures |27] and [29]). 
This suggests that the T-p relationship alone cannot account 
for the broader absorption lines close to galaxies. 

Another possible contributor to increasing bd near galax- 
ies is a change in the ionization at a given A^hi- If the galaxy 
itself produces a significant fraction of the local ionizing ra- 
diation field (i.e., in excess of that of the metagalactic UV 
background), then A^hi/A^h would be smaller near galaxies. 
The fraction of the ionizing photons which escape the ISM 
of the galaxy is small (Nestor etal. 2011, Steidel et al. in 
prep), but their number, compared to the UV background is 
likely significant within ^ 100 pkpc. In this case, the T-p 
relation could hold, but the mapping of A^hi to p would be al- 
tered. If this were the case, we would expect the locations 
of the asterisks in Figures |27] and |29] to be consistent with a 
leftward shift of the black bars. However, this scenario is at 
odds with the observation that the median bd of absorbers with 
10"* <A^Hi< lO"" cm-2 and 100 < Aran< 200 pkpc is higher 
than the median bd for any bin in A^hi observed in the full ab- 
sorber catalog (Figurel27b. The same is true of absorbers with 
D30 < 400 pkpc and 10'"* <A^hi< lO"" cm'^ (Figure|29ll. 

Assuming that overdensity is not the main cause of the in- 
creased bd, we consider other plausible processes which could 



26 



Rudie et al. 



80 




Dju < 200 


E 






if 60 
sz 






p 

5 40 
c3 

Q. 
Q. 






O 

Q 20 







13 



14 15 

log(NH 



16 



17 18 




13 14 15 16 17 IE 
log(NHi) [cm '] 




" 80 

03 

E 

60 



40 



Q 20 





D30 < 800 


* * 









13 14 15 16 17 18 

log(NH,) [cm T 



13 14 15 16 17 IE 

log(NHi) [cm '] 



Figure 29. Same as Figure|27]but now divided into bins of 3D distance, Z)3d, and witli cumulative bins in Dm. From Figure l28] one notes that the majority of 
the signal results from D}£,< 400 pkpc (shown in this Figure in the upper two panels). Notably, absorbers near galaxies have higher at nearly all values of Whi 
at these impact parameters. Because the signal from all bins is combined, higher foj at fixed A'hi is evident for absorbers with D}£,< 600 pkpc. 



lead to an increase in the temperature or turbulence of the 
gas. In principle it is possible to disentangle the source of the 
broadening (turbulent vs. thermal) by comparing the widths 
of absorption lines arising from ions with different atomic 
weights. The thermal broadening component, bd{T), depends 
on atomic mass, while the turbulent broadening component 
will remain constant so long as the ions reside in the same 
gas. While the HIRES QSO spectra in our sample include 
metal absorption lines which could be used for this purpose, 
the analysis of the metal lines is beyond the scope of this pa- 
per and will be presented elsewhere. At present, we consider 
the implications of both possible effects (temperature and tur- 
bulence). 

Referring to Figure |28] there is an increase in the median 
from 27 km s"' to 35 km s"' . Assuming purely thermal 
broadening, these values would indicate T = 4.4 x 10^ K and 
7.4 X 10"^ K, respectively, or an increase in temperature of 
68%. If we consider the absorbers with ^ 40 km s~' as in 
Figure |29l these suggest a temperature of ^ 10^ K. If instead 
the increase in line width is due to turbulence, an excess turbu- 
lent velocity of 20-30 km s"' would account for the departure 
from the bii= 27 km s"' global median. In either case, it is cu- 
rious that the elevated Doppler widths occur in gas with 100 
< Dtran< 300 pkpc but are not detected in the bin from 0-100 
pkpc. While the non-detection at the smallest impact param- 
eters is likely not significant, clearly, whichever mechanism 
is responsible for this increase in width is acting significantly 
outside of the typical galaxy virial radius (^ 90 pkpc). 

Unless AGN are involved (galaxies with detected AGN sig- 
natures were not used in our analysis), excess photo-heating 
for material closer to galaxies seems unlikely, since the UV 
radiation field from local stellar sources is likely to be (if 



anything) softer than the metagalactic background believed 
to dominate at an average IGM location. 

There are at least two plausible processes which could 
increase b^ in the CGM of galaxies: galactic winds (out- 
flows) and baryonic accretion. Accreting gas can be heated 
by shocks forming as the gas falls into galaxy halos. The 
galaxies in our sample have halo masses comparable to 
the theoretically expected transit ion mass between "cold 
mode" and "hot mode" accretion (Ivan de Voort et aLllMTTbt 
iFaucher-Gi guere et all [20 1 U iKeres et all 120091; lOcvirk et all 
2008); most of these authors predict that both hot and cold 
accretion occur in such halos at z > 2. The virial tempera- 
ture for galaxies with halos of mass ~ 10'^ Mq is expected 
tobeTvir 10^ K (Ivan de Voort et alJl201 Ibl) . and lKeres etaP 
(120091) suggest that this hot gas may fill the volume surround- 
ing galaxies to a few ryir- The H I absorbers in our sample are 
unlikely to trace gas this hot, however gas cooling from this 
temperature would be detectable. 

The interaction of cool dense gas with either a hot 
halo or a rarified wind fluid is also expected to produce 
turbulent boundary layers via bo th Rayleigh-Taylor^'^ and 
Kelvin Helmholtz^^ instabilities (iKeres & HernquistI 120091; 
Faucher-Giguere & Keres 201l|). These boundary layers are 
also where metallic ions such as O VI are expected to be most 
abundant. 

Galactic super winds can naturally explain either thermal 

Rayleigh-Taylor instabilities would form in the case of a lighter, less- 
dense fluid lying deeper in the gravitational potential of the galaxy than a 
denser cold-stream. 

Kelvin Helmholtz instabilities would form in the condition where the 
dense stream was moving with respect to the hot medium. This is especially 
plausible in the case of a fast moving galactic wind or a filament moving at 
the free-fall velocity. 
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or turbulent broadening, potentially to large distances away 
from galaxies. The favored mechanism for large scale winds 
is shocks from supernovae which could easily increase the gas 
temperature of the ambient medium. Although the distance 
to which galactic winds propagate is not well const rained, 
there i s evidence that the y expand to at least rvi, (Steide Tet alj 
I2010I) . IShen et aT] (1201 ll) recently considered high-resolution 
"zoom-in" simulations of a single Lyman Break Galaxy^* at 
z = 3. In their model, the LBG itself (as opposed to its satel- 
lite galaxies) enriches the IGM with metals from supernovae- 
driven winds to 3rvir (see their Figure 9). Kollmeier et al] 
(j2006) considered the effect of winds on the CGM of star- 
forming galaxies at z = 3 and found that winds could have 
a profound effect on the temperature of nearby gas: vary- 
ing the energy released by supernovae in their simulations 
by a factor of 5 changed the temperature of gas close to 
galaxies by ~ 250% (see the ir Figure 10). More recently, 
Ivan de Voort & Schavd (120 1 11) considered the physical prop- 
erties of gas surrounding galaxies with halo masses Mdm = 
lO'^ M0 at z = 2 as a function of the kinematics of the gas (i.e. 
inflowing vs. outflowing). They found a radial temperature 
profile of outflowing gas qualitatively consistent with the ob- 
served increase in b^. Outflows are also predicted to drive gas 
turbulence outside of virial halos through the conversion of ki- 
netic energy from supernova driven b ubbles into random mo - 
tions through gas instabilities (see e.g.. lEvoli & Ferraral201 Ih . 

While the exact cause of the enlarged line widths is not ap- 
parent, it is clear that the physical state of the gas within ^ 
300 pkpc is markedly different from that in "random" loca- 
tions in the IGM - suggesting that galaxies and their poten- 
tial wells deeply affect their surrounding CGM. This fact may 
introduce significant complications to studies of the general 
IGM which will be discussed further in i 



8. GENERAL DISCUSSION 

In ^-|2]we presented measurements and plausible inter- 
pretations of the properties of neutral hydrogen gas surround- 
ing star-forming galaxies at high z. Here, we speculate on the 
possible implications of the observed trends. 

8.1. Gaseous "Zones" around Galaxies 

The distribution of strong absorbers (Figure |9]l and of the 
incidence of absorbers (Figures [16] and [T7]i suggest three dis- 
tinct "zones" in the gaseous envelopes surrounding galaxies. 

1. The first zone corresponds to the volume within Z)tran< 
300 pkpc. For scale, the virial radius of a typical galaxy 
in our sample is ^ 90 pkpc. Inside 300 pkpc, the values 
of the max(A^Hi) statistic are significantly elevated and 
rise toward the position of the galaxy (Figure|9]l. Figure 
[16] shows rising Tyabs with decreasing Dnan for the three 
bins within 300 pkpc. Similarly, at 300 pkpc, ryabs for 
absorbers withA^Hi > 10'"*^ drops sharply. InFigurelTTl 
again at 300 pkpc one sees a rapid drop in the incidence 
of all absorbers with A^hi > 10'"*^ cm"^. 

2. The second zone lies between 300 pkpc < Dtran^S 2 
pMpc. In Figure [9] we observed that max(A^Hi), while 
high compared to the random distribution, plateaus at 
a roughly constant value within this zone. 77abs also 
plateaus over the same range as shown in Figures [T6] 
and[T7] 
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Figure 30. The fraction of all H 1 systems in our QSO sightlines that arise 
within I Av| < 300 km s"' and Dtran < 300 pkpc of the position of a galaxy 
in our spectroscopic galaxy sample (lighter red histogram), as a function of 
Whi ■ The darker (blue) histogram indicates the fraction after correcting for the 
fact that not all galaxies in the photometric sample have been observed spec- 
troscopically. These corrected points thus reflect the fraction of absorbers 
arising within the CGM of galaxies that meet our photometric selection cri- 
teria. 



3. The third zone corresponds to ftran^ 2 pMpc. The tran- 
sition between the second and third zone is not a sharp 
feature; rather between ^2-3 pMpc, the distribution 
of absorbers as seen in the max(A^Hi) and ?7abs statistics 
drops to become consistent with (or below that of) the 
IGM median. 

Three distinct zones are also evident in the velocity distri- 
bution of absorbers. 

1 . The first zone, corresponding to the strongest peak in 
the velocity distribution, is within |Av| < 300 km s~' 
(Figures [6] and |7]l . 

2. The second velocity zone includes gas to |Av| « 700 
km s"', which encompasses the full extent of the excess 
A^Hi (Figures [5] and O. 

3. The third zone is |Av| > 700 km s"', where the distri- 
bution with respect to galaxies is consistent with that of 
random locations (Figures|5]and|6]l. 

Note that in Figure [7] the median A^hi of absorbers with 
^'tran < 300 pkpc exhibits a sharp drop-off at |Av| > 300 
km s~' suggesting that absorbers at Dtran > 300 pkpc are likely 
responsible for the excess absorption at 300 < |Av| < 700 
km s~ 
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8.2. Defining the CGM 

We adopt a working definition of the "circumgalactic 
medium" (CGM) as the region within | Av| < 300 km s"' and 
^'tran < 300 pkpc of a galaxy. Figure [30] shows (in red) the 
fraction of absorbers from the full absorber catalog which fall 
within the CGM of a galaxy in our spectroscopic sample as a 
function of A^hi- 

Our spectroscopic galaxy sample is 70% complete with re- 
spect to the photometric parent sample within 300 pkpc of the 
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line-of-sight to the QSOs.^^ The blue histogram in Figure [30l 
has been corrected for this incompleteness assuming that un- 
observed galaxies have the same overall redshift distribution 
and similar CGM to those that have been observed. Figure[30l 
shows that ^ 20% of all absorbers with log(A^Hi) > 14.5 and 
^ 40% of those with log(A^Hi) > 15.5 arise in the CGM of 
galaxies that meet our selection criteria. Noting that the co- 
moving number density of identically selected galaxies (with 
the same limiting magnitude ofTl = 25.5) is $ = 3.7 x 10"^ 
cMpc-3 jReddv et al.l 120081) . the CGM of LBGs at (z) = 2.3 
contains only ~ 1.5% of the Universe's volume but accounts 
for nearly half of the total gas cross-section for A^hi > 10'^ 
cm"^. This is perhaps surprising given that our sample in- 
cludes only galaxies with Luv > 0.25Lyy and the faint-end 
slope of the UV luminosity function is very steep at z ^ 2.3; 
however, it is consistent with previous results based on strong 
metal-lin e absorption around similar galaxies at cornparab le 
redshifts (" Adelberger et"ani200l l2005i iSteidel et al.ll2010l) . 

The situation changes substantially for A^hi ~ lO''*^ cm"^; 
we have seen that these lower A^hi absorbers are more loosely 
associated with galaxies in our sample, though their incidence 
is enhanced by ~ 10-20% relative to average locations in the 
IGM (e.g. see the top panels of Figure [T6]l. 

Stated another way, we have found that low and high- col- 
umn density absorbers cluster differently with high-z, UV- 
bright galaxies. We showed (Figure [TSIl that absorbers with 
log(A^Hi) > 14.5 exhibit a much stronger peak in incidence 
towards the positions of galaxies than do the absorbers with 
lower column densities. Similarly, in maps of the incidence 
(Figures |2T| - 124] | low-column density absorbers show only 
large scale correlation with the positions of galaxies while 
high-column density absorbers again showed a strong peak 
near galaxies both along the line of sight and in Dtran. 

Collectively, these lines of evidence point to a circumgalac- 
tic zone defined by the boundaries of |Av|< 300 km s"' and 
Daan< 300 pkpc and dominated by absorbers with log(A^Hi) 
> 14.5. They also suggest that only systems with log(A^Hi) 
^ 14 cm~^ are confidently "IGM"- as A^hi increases beyond 
this limit, the likelihood that the gas lies within the CGM of a 
relatively bright galaxy increases very rapidly. 

8.3. Small Scale Distribution ofNni 

In an effort to understand the physics responsible for the 
distribution of A^hi as a function of transverse distance from 
galaxies, we compare the data to two theoretical models. First 
we consider the possibility that the H I column density traces 
the dark matter halo density profile. 

The typical dark matter halo hosting a gala xy in our sam- 
ple is theoretically expected to follow a NFW dNavarro et al.l 
I I997h density profile: 

, s Pciz) ^NFW 

p(r) = -2 (20) 

vir \ vir/ 

where (5nfw is a normalization defined as 
200 



ONFW ■ 



3 ln(l+c)-^ 



(21) 



and is the critical density, c is the concentration parameter 
of the dark matter halo which we take to be c = 4 following 

Our absorption line catalogue is essentially 100% complete for ab- 
sorbers with log(A'Hi) 13. 
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Figure 31. The max{Nn[, 300 km s"') statistic, converted into a univer- 
sal overdensity using Schaye ( 2001) plotted against transverse distance to 
a galaxy on a logarithmic scale. The same data shown in Figure [To] with 
two models over-plotted for comparison. Here we use 100 pkpc bins within 
1 pMpc and 400 pkpc bins at larger Dtran to simplify the plot. Shown in the 
(red) dashed Hne is the NFW radial density profile plotted as a universal over- 
density. See text for more details. Shown in the light (orange) dotted , light 
(green) dash-dotted, and dark (purple) dash-dotted lines are the average dark 
matter profiles from the MultiDark simulation for halos with Mdm > lO"^ 
M0. The light (orange) dotted curve is the median value (across the ha- 
los considered) of the average density in concentric spherical shells of radius 
Dtnu- The light (green) dash-dotted curve is the median (across the halos con- 
sidered) of the average dark matter density in concentric cylindrical shells of 
radius Dtran and length ±300 km s"'= ± 1200 pkpc. The dark (purple) dash- 
dotted curve is analogous to the green one, but the velocity window is ±700 
km s"' . 



iDuffv et al] (120081) . and rvii is the virial radius for halos of 
average mass lO'^ Mc?) (91 pkpc; Trainor & Steidel, in prep; 
iConrov et al1l200l lAdelberger et al.ll2005bl) . 
We define A, the matter overdensity, as: 



A=^ = 



p nM(i +zy Pciz = 0) 
Noting the definition for the critical density: 



Pciz) = 



(22) 



(23) 



we can then express the NFW profile as a function of over- 
density, rather than density: 



Anfw('") = 



1 



(5nfw 



hHz) 



nMil+zY 



vir V '^vir / 



(24) 



Using Jeans scale arguments. ISchavd (120011) defined a scal- 
ing function which relates the gas overdensity iph /p^) to an 
expected value of A^hi. While on the smallest scales where 
the growth is highly non-linear it is certainly not the case that 
Pb is equal to A, the matter overdensity (which is domi- 
nated by dark matter), the expectation is that they should mir- 
ror each other well at low overdensity (Awl- 10) where the 
growth is yet to become significantly non-linear. We there- 
fore assume /Oh/p^, ~ A, in order to compare our measured 
max(A^Hi) to theoretical density profiles. For the cosmology 
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Figure 32. Same as Figure |31|b ut using a different normalization of tlie 
photoionization rate. In Figure l3l] the data values were converted into A as- 
suming = 0.5 as suggested by Faucher-Giguere et af (2008). Here, the 
data values are converted into A assuming = 0.27 which bring them into 
the closest agreement with the MultiDark spherical density profile measure- 
ments. 

used in this p aper, at t he red shifts of our sample, the scaling 
function from lSchayd (120011) is: 
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where T4 is the gas temperature in units of 10"^ K and is the 
hydrogen photoionization rate in units of 10 "'^ s~' w h ere th e 
normalization is taken from Faucher-Gigue re et al.l ( 120081) . 
Using this relation, we can convert our max(A^Hi) statistic into 
A which can then be compared to the rescaled NFW profile. 

The values after converting max(A^Hi) to A are shown in 
Figure |3T] The (red) dashed curve is the NFW density profile 
(Equationl24li. The NFW halo has a radial distribution similar 
to that of the A inferred from our max(A^Hi) statistic for small 
^'tran, but falls Well below the Jeans-scale-inferred overdensity 
for Dtran ^ 300 pkpc. This is expected since the NFW pro- 
file is measured in simulations for dark matter (DM) within 
the virial radius. In the following section, we compare to the 
DM density on larger scales using DM profiles drawn directly 
from simulations. 

8.4. Large Scale Distribution ofNm 

Considering the distribution of gas around galaxies on Mpc 
scales, it is interesting to note the similarity between the scale 
over which circumgalactic gas has higher density than an av- 
erage location in the IGM (Dt,an S 2 pMpc) and the galaxy- 
galaxy autocorrelation scale length recently measured from 
the same galaxies, ro = (6.5 ±0.5)/;"' cMpc (~ 2.8 pMpc at 
(z) = 2.3; Trainor & Steidel, in prep). It would not be surpris- 
ing if they were closely related. 

Most of the recent work on the dark matter density profiles 
at several ryir and beyond have focused on the local universe. 
Since we are interested in comparing the gas-inferred density 
profiles with an average dark matter profile at z ~' 2-3, we 
created median DM density profiles using the MultiDark sim- 



ulation dKlvpin et al.ll201 iT) . The radial density profile within 
3 pMpc was computed from 100 halos with Mdm > 10" 
Mq^*^. Here we consider the dark matter density profile mea- 
sured in two ways presented in terms of universal overdensity 
as shown in Figure [3T| (1) First, we consider the average den- 
sity in concentric spherical shells of radius Dtran centered on 
the position of each halo. We take the median of this profile 
across the 100 halos considered (light orange dotted curve) 
and compare it to the max(A^Hi, 300km s"' ) statistic converted 
into overdensity using equation |25] (2) We also consider the 
average dark matter density profile computed within concen- 
tric cylindrical shells of radius Aran and length Av= ±300 
km s"' (±1200 pkpc) as shown in the light (green) dash- 
dotted curve. We consider the effect of the velocity window 
used as well - the dark (purple) dash-dotted curve is analogous 
to the green curve, but with a cylinder length of ±700 km s"' . 

The spherically-averaged radial density profile is well- 
matched by the NFW profile on small scales as expected (Fig- 
ure [JTT l. After 200-3 00 p kpc (the same scale as the first 
"zone" of the CGM, 38. Il l, the spherically averaged radial 
density profile flattens compared with NFW. From ^200-300 
pkpc onwards, the inferred density declines smoothly, falling 
to the mean density of the universe [log(A)=0] at Diaa> 3 
pMpc. 

The cylindrically averaged density profile significantly di- 
lutes the signal in the inner bins of Dtran suggesting that the 
max(A^Hi) statistic traces more of a radial distribution than a 
"line-of-sight" distribution on small scales. On larger scales, 
the cylindrical and spherically averaged density profiles look 
quite similar Between Dtran ~ 1.5-2 pMpc, the two cylin- 
drically averaged curves become consistent with each other 
This suggests that on scales of ^ 2 pMpc the density along a 
sightline between ±300 km s"' and ±300 to ±700 km s"' is 
similar - as one would expect far from galaxies as the density 
approaches the universal mean. 

The density inferred from the H I measurements is 
markedly higher on ^ pMpc scales than is predicted by ei- 
ther MultiDark DM density profile. There are three possi- 
ble causes for this discrepancy: (1) It may be the case that 
the max(A^Hi) statistic traces high-density substructure while 
the radial density profiles drawn from the simulations are (by 
construction) smoothed - averaging over these structures. The 
assumption of Jeans smoothing intrinsic to the arguments of 
ISchavd (120011) should minimize the discrepancy; however, it 
is likely still that substructure may affect our max statistics . 
(2) One or more of the assumptions used in the lSchav^ (1200 lb 
conversion between A^ei and A could be invalid. For exam- 
ple one of the more uncertain quantities required in the con- 
version is the H I photoionization rate, F.^' In Figure [32] 
we consider the change in the value of V needed to bring 
the measured A[max(A^Hi, 300km s"')] into agreement with 
the spherically averaged dark matter density profile. This 
plot suggests a normalization ~40% lower than the estimates 
from lFaucher-Giguere et alj (120081) ."^° (3) Another possibility 

The halo mass was inferred from the number of particles found to "be- 
long" to a halo using the Friends-of-Friends algorithm. The halo-mass cut 
of 10" * M0 is the same mass threshold which reproduces the clustering of 
galaxies in our spectroscopic sample, as described in Trainor & Steidel (in 
prep). The mean and median halo masses of galaxies with this halo-mass 
threshold are 10'^'^ Mq and 10'^ Mq, respectively. 

refers to the photoionization rate in units of 10"'^ s"' . 
4" The temperature of the gas is another quantity that could be adjusted. 
The temperature assumed currently is lO'* K. To match our observations we 
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is that Jeans scale arguments do not hold in detail on these 
scales."^' 

8.5. The Implications of the Distinct Physics of the CGM 

The final point is a note of caution regarding the effect 
the CGM of luminous galaxies may have on measurements 
which rely on "average" locations in the IGM. In addition 
to the rising A^hi and 77abs observed in the CGM of galaxies, 
we have shown that the physical state of the absorbing gas in 
this region is also distinct from the general IGM. The Doppler 
widths, bi, of individual absorbers are significantly larger in 
the CGM zone, even when compared fixed Nhi- Further, 
absorbers with A^hi> lO''*'^ cluster tightly with galaxies (Fig- 
ures [16] [T7] and l22l i. and these absorbers are found in the 
CGM zone with very high frequency (Figure [30]l. 

Taken together, these observations suggest that the physics 
of the "IGM" which has long been leveraged for measure- 
ments varying from the metagalactic background to the mat- 
ter power spectrum, may be significantly affected by the bary- 
onic physics of galaxy formation; particularly, low-resolution 
IGM studies dominated by H I absorbers ^saturated in Lya 
(A^Hi> 10'"*^ cm~^) are likely to be affected by the presence of 
a nearby galaxy. Further observational and theoretical exam- 
ination of this subject is crucial to a full understanding of the 
measurement power of the IGM. 

In addition, studies of the distribution of metal absorption 
lines in the "IGM" often consider metallic species associ- 
ated with HI absorbers of log(A^Hi) > 14.5 (Co wie et al. 1995.: 
Spngaila & Cowie 1996; Ellison et al. 2000t ICarswellet al J 
2002; .Simcoei,201 1) . To a large degree, this is a practical 
choice as the associated metal absorbers occurring in sub- 
solar metallicity gas must be found and measured in spec- 
tra of finite S/N. However, the properties of the CGM dis- 
cussed above suggest that, as many authors have speculated, 
these measurements more directly probe the metallicity of the 
CGM than that of the IGM. Observations at z > 5 generally 
cannot measure the A^hi associated with metal systems due to 
the line density of the forest, but they are likely to suffer from 
the same effect, adding credence to the idea that the recent 
discovery of the high-z do wnturn of C IV (.Rvan-Weber et all 
Il006, 2009; Becker et alJ l2009l IMThI: LSimcoe et all 120 1 Ih 
likely also traces CGM gas. 

The observations of circumgalactic and intergalactic gas 
and its association with forming galaxies have now progressed 
to the point that more stringent comparison to theoretical ex- 
pectations is possible. Even using the comparatively "blunt" 
tool of gas covering fraction, we have shown 05.21 ) that the 
distribution of HI gas around galaxies exceeds the expec- 
tations of the models to which they have been compared, 
with the discrepancy increasing as A^hi decreases below 10'^ 
cm"^. It is clear that a combination of high-resolution (to 
capture the small-scale baryon physics) and cosmological (to 
provide sufficient statistical leverage) simulations will be re- 
quired to make use of the more detailed physical measure- 
ments emerging from the observations: e.g., gas-phase kine- 
matics, multi-phase velocity/density structure, redshift-space 

would need to significantly lower the temperature assumed to T < 10^^ K 
which is unlikely to be the case. If anything, 10* K is likely a lower limit on 
the temperature of the IGM at z ~ 2 (see e.g., Becker et al. 2011a). 

We tested the agreement between the aforementioned models and our 
data converted usin g the Fluctuating Gunn-Peterson Approxima tion (FGPA; 
IRauch et alj 119971) . finding that the disagreement is larger. IRakic et al] 
<2011bl) also compare several methods for converting between H I optical 
depth (akin to A'hi) and overdensity (see their §4.5). 



distortions, gas temperatures, turbulence, the distribution and 
mixing of metals, all as a function of galaxy properties. In 
order to match the real universe, simulations will need to re- 
alistically model the physics of feedback from star formation, 
supernovae, and AGN activity and its effect on gas-phase ac- 
cretion from the CGM. 

It will be important for future simulations to have sufficient 
dynamic range to make predictions about the distribution of 
neutral hydrogen with 13 < log(A^Hi) < 20. The physical pic- 
ture advocated by most modern theoretical treatments of the 
"IGM" is that the Lya forest traces gas of modest overdensity 
within the filaments of the "cosmic web". The high cover- 
ing fractions found in this study (e.g., the relatively large fc 
of saturated HI extending >2 pMpc from ^ 10'^ Mq dark 
matter halos; © suggest that if this picture holds, the size of 
these filaments may remain quite large even as they approach 
the positions of high-bias halos. These initial comparisons 
suggest that the topology of the neutral hydrogen gas may be 
less "filamentary" (at all column densities) than current sim- 
ulations of the "cosmic web" and "cold accretion" predict. If 
such disagreement between simulations and observations are 
common and persist they may become key to a deeper un- 
derstanding of the baryon physics that ultimately controls the 
process of galaxy formation. 

9. SUMMARY 

In this paper we have presented a detailed study of neu- 
tral hydrogen in the IGM surrounding 886 high-z star-forming 
galaxies using spectra of background QSOs from the Keck 
Baryonic Structure Survey. We draw conclusions from the 
analysis of 15 sightlines to hyper-luminous high-z QSOs us- 
ing high-resolution (7 km s"'), high-signal-to-noise ratio (50- 
200) spectra. This study constitutes the largest absorption 
line catalog ever created at these redshifts through a com- 
plete Voigt profile decomposition of the Lya forest for all 15 
sightlines. The high S/N and excellent UV/blue-wavelength 
coverage of the HIRES QSO spectra allow precise determina- 
tion of the column densities of saturated H I absorption sys- 
tems whose positions we have found to correlate strongly with 
those of the galaxies in our sample. This sample is also con- 
sidered by Rakic et al. (2011b) who analyze the pixel statis- 
tics of the H I optical depth. 

In this work, we study the absorption patterns of H I gas 
in the velocity and spatial locations surrounding each galaxy 
redshift. We have examined the correlations among column 
density (A^hi), covering fraction (/c), incidence (jyabs), absorp- 
tion line width (ba), galactocentric transverse distance (Duan), 
velocity offset (Av), and 3D distance (D^u)- Our principal 
results are as follows: 

1. The redshifts of H I absorption systems in the QSO 
spec tra c orrelate strongly with the redshifts of galax- 
ies 04.1I ). The distribution of relative velocities ( Av) is 
peaked at the systemic redshift of galaxies in our sam- 
ple (Figure |5]). The full excess absorption near galax- 
ies falls within ±700 km s~' with the majority of the 
excess A^hi (i-e. the higher column density absorbers) 
falling within ±300 km s"' (Figures |6] and [T). 

2. The column density (A^hi) of absorbers is strongly cor- 
related with the impact parameter to a galaxy (Dtran, 
94.21 ). We consider the highest-A^ni absorber within 
±300 km s"' [max(A^Hi)] of each of the 10 galaxies 
within 100 physical kpc and find that the median value 
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is 3 orders of magnitude higher than at random places in 
the IGM (FigurefTOli. A trend of rising column densities 
as a function of decreasing Duan is observed within the 
inner 300 physical kpc (1 cMpc) surrounding galaxies 
(Figures 121 [lOl and[Il). 

3. The column densities (A^hi) of absorbers as measured 
through the max(A^Hi) statistic plateau in the range 300 
pkpc < Aian< 2 pMpc (Figures |9] and [ini EZB. In- 
terestingly, the plateau is elevated in A^hi with respect to 
random places in the IGM by 0.3 - 0.6 dex. For AranJ^ 
2 pMpc, the values of max(A^Hi) decline and become 
consistent with random places in the IGM. 

4. The median A^hi as a function of 3D distance"^^ 
smoothly declines with increasing D^u and becomes 
consistent with random places in the IGM at D^u ~ 3 
pMpc (Figure [H ^43[ . 

5. The incidence"*^ (??abs) of absorbers of all A^hi is higher 
within Dtian ^ 2 physical Mpc of the position of a 
galaxy than at random places in the IGM (Figure [161 

For DtianJ^ 2 pMpc, 77abs is consistent with that of 
the random IGM. 

6. For absorbers with log(A^Hi) > 14.5 cm"^, there is a 
strong peak in the incidence, 77abs, close to galaxies at 
Dtran < 300 pkpc (Figure[T6] No such peak is seen 
in the distribution of lower-column density absorbers. 
This suggests that absorbers with log(A^Hi) > 14.5 cm"^ 
are more directly related to galaxies than absorbers of 
lower A^Hi- 

7. The incidence, ?7abs, of absorbers with A^hi > lO'"^'^ 
cm"^ and | Av|< 300 km s~' shows rising values within 
£>tran < 300pkpc. From 300 pkpc < Dtian< 2 pMpc, the 
ryabs roughly plateaus with ryabs ^ 0.5. After ~2 pMpc, 
the distribution becomes consistent with random places 
in the IGM (Figure[T7l EB. 

8. Within 300 pkpc and ±350 km s"' of a galaxy, the prob- 
ability of intersecting an absorber with log(A^Hi) > 14.5 
is > 4 times higher than at a random place in the IGM. 
(Figure[T8]and Tableg] EB. 

9. Maps of the incidence, r^abs, of low and high-A^Hi ab- 
sorbers as a function of Duan and Dnubbie allow for the 
measurement of redshift anisotropics (® . In Iow-A^hi 
gas [log(A^Hi) < 14.5], on large scales (pMpc) the distri- 
bution along the line of sight is compressed compared 
to the distribution in the transverse dire ction. This is 
likely the signature of gas infall e.g. the Kaiser' (Il987h 
effect (Figure |24l). In the high-A^Hi gas [log(A^Hi) > 
14.5], we find the "finger of God" effect, the elongation 
of the distribution along the line of sight compared to 
the transverse distribution (Figure l22l). Taking into ac- 
count the effect of our redshift errors, the data are con- 
sistent with peculiar velocities of ^ 250 km s"'. These 
conclusions are in agreement with those presented in 

The 3D distance, Djq, is computed using tlie quadratic sum of tlie pliys- 
ical impact parameter between tlie galaxy and tlie line of sight to the QSO 
(Aran) and the line-of-sight distance one would calculate assuming the Ai' of 
each absorber was due only to the Hubble flow, ^Hubble- 

"^^ The number of absorbers per galaxy with | Ai'|< 300 km s"' and in a 
given bin of Dtran- 



iRakic et aTl (1201 Ibh who used pixel optical depth anal- 
ysis to demonstrate that the compression along the line 
of sight on large scales is highly significant (> 3(7). 

10. The Doppler widths (bi) of individual absorption lines 
in the HIRES spectra rise as the galactocentric distance 
is decreased (©. The majority of the increased width 
is found within |Ay|< 150 km s"' and 100 < Aran< 
300 pkpc (Figures l26l and |28]). The larger line widths 
could be caused by an increase in the gas temperature, 
the turbulence of the gas, or both. The median bd close 
to galaxies are consistent with a gas temperature ^10^ 
K, or a turbulent contribution of ^20-30 km s"' . 

11. Based on these measurements, we suggest a working 
definition of the circumgalactic medium (CGM) to be 
all locations within 300 pkpc and ±300 km s"' of a 
galaxy. We find that the CGM of our completeness cor- 
rected galaxy sample can account for nearly half of all 
absorbers with log(NHl) > 15.5 (FigurepOl ^211. No- 
tably, the CGM of these galaxies comprises only 1.5% 
of the universe's volume at these redshifts. 

These findings not only define relevant scales of the CGM 
of luminous galaxies at redshift z ^ 2-3, but also demonstrate 
that much of the moderately high-A^ni portion of the IGM 
likely originates in the regions surrounding luminous galax- 
ies where the baryonic physics of galaxy formation appear to 
affect the physical state of the gas. This point may have im- 
portant consequences for past and future studies of the "IGM" 
and should be considered in greater detail. 

Our understanding of the physical properties of gas sur- 
rounding star-forming galaxies will greatly benefit from anal- 
ysis of metallic absorption features in the HIRES spectra. 
These absorption lines probe the ionization state and metal- 
licity of the gas and will provide direct constraints on the total 
mass in hydrogen and in metals that surround star-forming 
systems. This work is underway and will be published in the 
second paper in this series. 

With this further analysis, we will have for the first time 
a quantitative snapshot of the baryonic interplay between the 
intergalactic medium and forming galaxies which is evidently 
central to the processes of galaxy formation and evolution. 
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